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Abstract 

Since the recent dissertation by Steffen Winter, for certain self-similar sets F the 
growth behaviour of the Minkowski functionals of the parallel sets := {x G 
: d{x, F) < e} as e t 0 is known, leading to the notion of fractal curvatures 
Cl{F), k G {0, The dependence of the growth behaviour on the fractal 

dimension s = dimF is exploited, and estimators for s and Cl{F) are derived. The 
performance of these estimators is tested on binary images of self-similar sets. 
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Introduction 


At the end of the 19th century, the common belief that every subset A C of 
euclidean space could be assigned an integer number called “the dimension” of A 
could not be supported anymore: Examples like the Cantor-set (termed after the 
German mathematician Georg Gantor) showed that geometric measure theory needed 
to be generalized, which was then done essentially by Felix Hausdorff. The theory 
around these abstract sets, however, did not affect the natural sciences until some 
50 years later, when strange attractors were discovered. Attractors are subsets of 
with the property that there is a dynamical system in which trajectories of particles 
converge to this set, and they are termed “strange” if their dimension is non-integer 
(fractal). 

This development continued with a general reconsideration of how accurately length 
and area of concrete objects can be determined; a good example is the question “How 
long is the coast of Britain?” (B. Mandelbrot). A commonly expected answer would 
be a 4 digit number of miles, which corresponds to the length of a rectihable curve 
that represents the coastline. Recall that measuring means adding up the lengths of 
the steps a divider has to take to run through the curve. The reciprocal value of the 
steplenghts is usually termed “accuracy”. As the accuracy tends to oo, the measured 
length will stay hnite for rectihable curves, and approach oo for non-rectihable curves. 
Now the coast of Britain cannot ever be measured with arbitrarily high accuracy, but 
extrapolation to inhnity suggests that the measured lengths might rather tend to 
oo than to a hnite valu^ An answer to whether or not this is right is impossible, 
of course, since the concrete coastline is something intrinsically diherent than an 
abstract mathematical curve, be it rectihable or not. 

Note that this example can easily be taken into two dimensions, too: Gonsider a 

^ To visualize this, compare the distances i) an airplane has to travel to view all of the coastline 
of Britain ii) a human has to walk to circle Britain with his feet touching the water and iii) an ant 
has to travel if it is to cross every pebble of Britain’s shore which touches the ocean (frozen in time). 
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1 cm^ square piece of smooth glass. Under an electron microscope a very irregular 
shape is revealed, the measurement of the area of which on this scale is yielding a 
value far greater than 1 cm^. 

The observations above all have in common that measurements of areas and lengths 
grow as the scale shrinks. Fractal (not necessarily integer-valued) dimensions turn 
out to describe the rate of this growth; conversely, the knowledge of the growth-rate 
translates into knowledge of the fractal dimension of a set. To the best of the author’s 
knowledge, so far the fractal character of subsets of has only been described based 
on (the growth behaviour of) measurements of one physical quantity only, i.e. either 
of volume, area or of length exclusively. However, for a big class of subsets A of 
one has the four important characteristics termed “Minkowski-functionals”: These 
are denoted by C^{A), C 2 (H), Ci{A) and Co(H) and have the physical interpretation 
volume, surfaee area, mean breadth and Euler-charaeteristie (i.e. the number of 
connected components minus the number of holes) respectively. The reason why, so 
far, fractal properties have been analysed on the basis of C(i{A) only, must be the 
following: Because of the complexity of fractals in dimensions that are higher than 
one, nothing was known about the behaviour of the characteristics Ck{A) for k < d—1 
as the scale shrinks, until very recently S.Winter made the hrst step into this direction: 
For a subclass of the class of self-similar sets, which are probably the best-understood 
fractals, he thoroughly described the limiting behaviour of these characteristics. Now 
that this behaviour is known for a (yet very small) class of fractals, the theoretical 
foundation is laid for the analysis of fractals, also based on Cd-i{A ),..., Cq{A). Thus 
the question that this thesis is devoted to can now be formulated, here in the most 
abstract way possible: 

If, on a shrinking scale, the simultaneous growth behaviour 
of all Minkowski-functionals is known, how can this additional 
information be translated into further knowledge about geomet¬ 
rical properties of the fractal set? 

Of course, one can only hope to hnd partial answers to this question, in a very 
restricted setting. We are only going to consider the class of sets which are subsets 
of self-similar, satisfy the open set condition and have polyconvex parallel sets. 
By the r-parallel set of A we mean the set Ar of all points whose distance to H is < r. 
The value r reflects the accuracy of a measurement of A: If a set H G existed in 
nature, (it would most probably be invisible and) measurements of the locations of its 
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points would have an incertainty of ±r in each direction, so real-world measurements 
of the Minkowski-functionals Ck{A) can be compared to those of Ck{Ar). As the 
accuracy tends to inhnity, i.e. as r 0, Winter has shown that Q^sets would exhibit 
the growth behaviour Ck{Ar) ~ Clu~^ (where s = dim A), and that there is some 
analogy of the prefactors Cl to the classical curvatures. The scope of this thesis 
will be to estimate the values Cl and s from given data Ck{Ar). The construction 
of A, the dilation of A by r and the measurement of Ck{Ar) were all performed by 
computer simulations, where the author strongly benehted from algorithms that were 
readily implemented in the GeoStoch ( |Dep07| ) library of the University of Ulm. 

This thesis is organized as follows: Part I describes basic facts of fractal geometry in 
chapter [T] and excerpts of the theory behind fractal curvature in chapter which is 
mostly a summarization of theorems taken from Winter’s dissertation |Winn5] that 
are necessary for chapter]^ In part II, chaptera short introduction to the problem 
of squeezing fractals into binary images is given, followed by chapter]^ which contains 
a review of methods that analyse fractal images. The author has found a great 
amount of literature on estimating fractal dimension, but nothing on estimates of 
the minkowski content, which would have been interesting with respect to chapter]^ 
Literature on the estimation of the dimension of graphs of functions is vast but has 
been ignored because it hts more in the scope of time series than image analysis. The 
main chapter, chapter introduces a joint estimator of fractal dimension and of the 
three fractal curvature measures for sets in Finally, in chapter]^ the performance 
of the estimator is tested on several sample images. 
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Fractal Geometry 
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Chapter 1 


Basics of Fractal Geometry 


In this chapter, notions from fractal geometry that are important for later use shall 
be summarized, and important theorems are cited. Proofs can be found in the books 
of Falconer |Fal9n] , Mattila |Mat95j and Rogers |Rog99| . 


1.1 Hausdorff-measure and Fractal Dimension of 

Sets 

The maybe most elementary characteristic that is classifying a geometric object is 
its dimension. In classical geometry, the topological dimension was most widely used: 
Wherever there exists a local isomorphism (a bijective mapping, both directions being 
continuous) from an open subset of {k G No) to the object, the object is said to 
have dimension k. 

In 1918, Hausdorff generalized this notion of dimension in such a way that a dimen¬ 
sion could be assigned to every Borel-set F G while the following demands were 
still met: 

• Monotonicity: E C F ^ dim(F) < dim(F) 

• Open sets: If F is a non-empty open subset of then dim(F) = d. 

• Geometric invariance: If / is a (non-singular) affine transformation, then 
dim(/(F)) = dim(F). 


7 








CHAPTER 1. BASICS OF FRACTAL GEOMETRY 


• Stability under countable unions: dim = sup^gj^ dim(ii^j). 

A set whose Hausdorff-dimension is non-integer or different from its topological di¬ 
mension is commonly termed a fractal, though no commonly used mathematical def¬ 
inition of the word “fractal” exists. 

We are now going to give a short introduction to Hausdorff-measures, which are the 
foundation of the Hausdorff-dimension, and we only consider the special case of the 
metric space X = with the euclidean metric d{-, ■). First we introduce a class of 
dimension functions: 

El := {h : [0, cx)] —)■ [0, ex) : h is right-continuous and non-decreasing, h{t) = 0 -v^ t = 0} 

Now for 5 > 0, A C X, h G El and |A| = sup{d{x,y),x,y G A} being the diameter of 
A, let 

{ oo oo 

J];i(|Ci|):Ac|JCj,|Ci|<d (1.1) 

i=l i=l J 

where the infimum is taken over all countable coverings {Q} of A with arbitrary 
subsets Cj C X. is not necessarily a Borel-measure, but 

n\A) := limn^(A) 

( 5 — 1-0 

is. It is given the name Hausdorff-measure with dimension function h. If 0 < 'H^(A) < 
oo then h is called the exact dimension function of A. 

Consider now the subclass 

EIq = {h{t) : h{t) = t°‘,a E (0, x)} 

of El. If h{t) = G EIo, it is common to put HA := H^ and to call it the a-dimensional 
Hausdorff-measure. The reason for this is that for positive integers d, H'^ and the 
d-dimensional Lebesgue measure \d are multiples of each other. 

For a < /3 we have H^{A) < S^~°'Hf{A), which yields the implications 

H“(A) < X ^ H"+%A) = 0, 

H“(A) > 0 ^ H^-%A) = X 

whenever 0 < e < a. An immediate conseqnence is that given a snbset A C there 
is one distinct valne s snch that 

f x for a < s 

0 for a > s 
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We write dimj:/(74) = s and say that A has Hausdorff-dimension s. Note that the 
three cases FL^^A) = 0, 0 < Fi^^A) < oo, FL^{A) = oo are possible. The second case 
is of special interest, and we shall call A C an s-set if 0 < FL^{A) < oo, i.e. if 
h{t) = F is an exact dimension fnnction of A. In a later chapter we will see that 
self-similar sets satisfying the open set condition (OSC) are snch sets. 


There are examples of sets A for which F is not an exact dimension fnnction for any 
s > 0, as almost all Brownian paths Et- In fact, dim/f(ST) = 2 but FC{Et) = 0 
almost surely; an exact dimension function would be e.g. t^loglog(y). 


The Hausdorff dimension is very useful in a mathematical framework, but it is not 


easy to calculate directly: The coverings from equation |1.1| consist of possibly in- 
hnitely many sets of varying diameter. Experimentalists prefer the following notion 
of dimension, as it is more easily calculated: 


Definition 1.1.1 For 5 > 0 and a hounded set A, let Ns{A) be the smallest number 
of sets of diameter not greater than 6 whose union covers A. The numbers 


dims (A) 


liminf and dims(v4) := limsup , 

<5^0+ -logd 5 ^ 0 + -log 5 


are called lower and upper Box-counting-dimensions (or lower and upper Minkowski- 
dimensions) of A. If they coincide, their common value dims(A) is simply called the 
Box-counting dimension (or Minkowski-dimension) of A. 


Box-counting dimensions satisfy dims (A) > dims (A) > dims^(y4), and unfortunately 
they are not stable under countable unions. However, in many cases as e.g. for 
self-similar sets satisfying the open-set-condition (OSC) (see Section 1.3), the box¬ 
counting dimension exists and is equal to the Hausdorff-dimension. 

Note that there is also the following representation of the above limits: 


Lemma 1.1.2 


dim s(T) = sup {t G M : Ns{F)6^ -P- co (5 — )■ 0)} (1.2) 

^b(F) = inf {t G M : Ns(F)d* 0 (6 ^ 0)} (1.3) 


In the limiting case s = t, limsio Ns{F)S^ can be either 0, a positive hnite value or 


oo. 
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Later (see Section 4.3) we will use the following alternative formula for the box¬ 
counting dimension, which goes back to Minkowski (and thus explains why there are 
two different names for the same dimension): 


Theorem 1.1.3 Let A (Z be a Borel-set, write for the Lebesgue-measure on 
and denote 

Ar := {x e : d{x,A) < r} 


the r-parallel set of A. Then 



dims (A) 

r 4 .o log r 

(1.4) 

dims (A) 

\ogC^{Ar) 

= d — hm sup — - - 

r.io logr 

(1.5) 


Proof: In the dehnition of the Box-counting dimension, Nr {A) can be substituted 
by Pr{A), the greatest number of disjoint r-balls with centers in A (see e.g. [Mat95| 
Ch.5]). Then the Lebesgue-volume of the parallel sets can be sandwiched via 

Pr{A)a{dy < C‘^{A{r)) < Nr{A)a{d){2rf, (1.6) 

{a{d) denotes the volume of the d-dimensional unit ball). After taking logarithms, 
divide by logr and consider the superior and inferior limits. □ 

Definition 1.1.4 Let A C such that s := dims(A) exists, and let the parallel 

sets Ar be defined as above. Then the limit 

M(A) := limC-^C^{Ar) 

rlO 

is called Minkowski-content if it exists. 

A result that will be used later is that for self-similar sets F satisfying the Open-Set- 
Condition we have 

0 < M(F) < oo. 

Dimitris Gatzouras proved this in |Gat00j (without the assumption of F having poly- 
convex parallel sets). 


^Here d(x, A) := inf(d(a:, o) : a S A) is the distance of the point x to the set A. 
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1.2 Dimensions of Measures 

Borel-measures can be viewed as a generalization of Borel-sets, in the sense that a set 
can only describe the shape of an object, whereas a measnre can describe its shape 
and its mass distribntion. Accordingly, measnres admit a richer “fractal theory” 
than sets. Here, however, only the framework for nnderstanding the Local Dimension 
Method in section 14.21 shall be laid. 

Let /i be a locally hnite Borel-measnre on and let B{x,r) denote the ball in 
aronnd x of radins r. 


Definition 1.2.1 The limit 

D'fx{x) := limsnp 

r^o+ (2r) 

is called upper s-density of p, at x. Similarly, the lower s-density of p at x is defined 
as 

r^s f \ 1 - • rTi.B{x,r)) 

D pix) := hmmf -—r-. 

— ^ r^o+ (2r)* 


The dimensions of a measnre are now dehned as follows: 


Definition 1.2.2 Let p be a Borel-measure on 

?’-5>o+ log r 

log r 

are called the lower and upper local dimension of p at x E 

These limits express the power law behavionr of p{B{x,r)) as r —0. If both lower 
and npper local dimension agree at x then p is called simple at x, and one writes 
dim p{x) for the common valne and calls it the local dimension of p at x. Note that 
for X G ontside the snpport of p one has dim/i(x) = oo, and that dimp(x) = 0 if 
X is an atom of p, i.e. if p{{x}) > 0. 

There is the following connection between the densities and dimensions of Borel- 


measnres: 
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Lemma 1.2.3 Let fi be a Borel-measure. 

dim/i(a;) = sup{t > 0 : D^^{x) = 0} = inf{t > 0 : D^fi{x) = cxd} 
dim/i(a;) = sup{t > 0 : D^fi{x) = 0} = inf{t > 0 : D^fi{x) = cxd} 


In other words, as a function of t, the lower density n{x) jumps from 0 to oo at the 
value dim/i(a;), whereas the upper density performs this jump at the lower dimension. 


1.3 Self-similar sets 


Self-similar sets are the best-understood fractal sets and serve as the starting point 
for the development of the theory of fractal curvature (see Chapter]^. 

A map / : X —X is called a contraction if 

T- / d{f{x)J{y)) 

Lip/ := sup- - -^- < 1, 

x^y d{x,y) 

and it is a similarity if for some constant r > 0 we have d{f{x),f{y)) = r d{x,y) 
for all x,y E X. An Iterated Function System (IFS) is a system S = of 

maps that are both contractions and similarities. 

Instead of (X, d) we are now going to consider the metric space (K'^,dH)- 

Definition 1.3.1 Let 

:= {K C : X is compact}. 


and let Ks be the 5-parallel set of K from definition 1.1.1 Then 

^ [0,cx)), 

(X, L) ^ inf{d > 0 : X^ D L and D X}. 
is called the Hausdorff-metric on W^. 


Given the IFS 5, consider the map 


N 

S:K^K, K ^\JSi{K). 

i=l 

Since the Si are Lipschitz-continuous, they map compact sets to compact sets, and 
since the union is hnite, S maps compact sets to compact sets, so S is well-dehned. 
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Lemma 1.3.2 S is a contraction on (K^dn) with Lip S = maxi<j<^r Lip^j. 


Proof: If (S'i(A))^ contains Si{B) for all i, then (Uilicontains 
and the same is true if A and B are swapped. Thus 


N 


N 




\i=l 


2 = 1 


< max dn iSi{A),Si{B)) < 

l<i<N 


max LipS'i duiA^B). 

l<i<N 


Theorem 1.3.3 For the Iterated Function System S let S : —?• be defined as 

above. Then there is a unique non-empty compact set F called the attractor of 
S that satisfies 

S{F) = F. 

Furthermore, if E & and E D S{E), then 

OO 

F = Pi S^{E). 

k=0 


Proof: This is a simple consequence of lemma 1.3.2 and Banachs Fixed-Point theorem 
on endowed with the Hausdorff-metric. 


For later use we shall distinguish two types of attractors / self-similar sets: 


Definition 1.3.4 The Iterated Function System S = {Si,..., Sn} with similarity 
ratios ri and the corresponding self-similar attractor F are called h-arithmetic z/h G M 
is the greatest number such that logr* G hZ Vz G {1, ... ,N}. If no such number h 
exists, S and F are called non-arithmetic. 


The following is a regularity property of iterated function systems: 


Definition 1.3.5 (Open Set Condition) An iterated function system S is said to 
satisfy the open set condition (OSC) if there exists an open set O such that 

SfiO) C OVz and SfiO) n S'j(O) = {(1 fori ^ j. 


All sets from section 6.1 satisfy the OSC. Now we can state the main theorem on the 
dimension of self-similar sets (see e.g. [MQ], theorem 9.3). 
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Theorem 1.3.6 Let S = {Si,..., Sjv} be an iterated function system with corre¬ 
sponding similarity ratios r*, and assume that the OSC holds. Furthermore, let s he 
the solution of = 1- Then the attractor F satisfies 

1. 0< n%F) < oo 

2. dimF = dim^F = dinisF = s 


Note that the OSC is guarantees that FL^^F) is bounded away from 0 in the above the¬ 
orem. There is a corresponding result for self-similar measures, also due to Hutchinson 
r |Hnt81] h 


Theorem 1.3.7 Assume the same conditions as in theorem Then there is a 

unique Borel-measure /r with = 1, called the S-invariant measure, such that 


and it satisfies 


M-) 


N 




a(-) 


n^{F) 


Corollary 1.3.8 The measure n is spreading the unit mass uniformly over all of F. 


(Here we use the word “uniformly” because /i is the restriction of the rotationally 
and translationally invariant s-dimensional Hausdorff-measure.) 

Proof: It suffices to show that supp /i = F. F is compact, so closed. Let x G F O [/ 
for some open set U. Then there exists a finite sequence ii, * 2 , ■ ■ ■ , 4 in {1, ..., N} 
such that := o ... o D F, so /i([/) > > 0. □ 

The specialty of this measure is that it spreads the unit mass uniformly on its support, 
which is a self-similar set of dimension s. By sampling points according to this 
probability distribution we can generate attractors of any iterated function system, 
and in fact this method has been used for the generation of the images from section 


via a non-uniform distribution can require many more points until every pixel of the 
attractor is finally covered. 

Before ending this section, the following result shall be noted, which justifies the 
Local Dimension method of section 14.21 


6.1 on a computer. The uniformity of the distribution is important, as sampling 
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Proposition 1.3.9 Let F be the attractor of the IFS S which satisfies the OSC. 


Then the self-similar measure /x from theorem 1.3.1 satisfies 


dim /i(a:) = s 


for all X G F. 


1.4 Generating Attractors of IFSs 


For later use (section 6.1) we now describe how to generate a binary image of the 
attractor F of an iterated function system. Chapter decribes the problem of dis¬ 
cretization more thoroughly, but for now it is sufficient to know that we represent a 
binary image of an attractor F C by all the pixels that are intersecting it (these 
pixels shall be called the pixels of F). 


Fix xo G and consider the sequence of points 


Xk-.= Sifixk-i), A; = 1,2,... (1.7) 

where R is an iic0 sequence of random variables with values in the set {1,..., 
having the probability distribution 


Pih =^)= rl 

(Recall that N is the number of similarities in the IFS.) In a few steps we will see 
that as fc —)■ cx), Xk runs through (almost) all the pixels of F. 


Lemma 1.4.1 Let c := max{rj, i G {1,..., . 

1. If d{xo, F) < N for some N > 0 , then d{xk, F) < c^N. 

2. If Xo G F, then Xk ^ F for all fc G N. 


Proof: Let F = {x G : d{x,F) < N}. Then Xk G S^{E), so d{xk,F) < 
dniS^iE), F) which is not greater than c’^dniE, F) = c^N by lemma 1.3.2 
shows 1. S{F) = F shows 2. □ 


This 


^independent and identically distributed 
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This means that if Xk is not already in F, at least it is approaching F exponentially 
quickly. If the distance is measured in pixelwidths, then practically this means that 
if the diameter |F| of F equals 3000, if xq is not further away from F than 30000 and 
if the similarities all have ratio not higher than 0.8, then xioo will be closer to F than 
10“® (pixelwidths). This should be enough to assert that almost no pixels outside F 
are marked black. Also, due to the contractive nature of the algorithm, the iterations 
of Xk are computationally stable, as calculation errors decay exponentially with k. 

Now consider the Borel-measure 

v{B) := lim —#{fc : k < n,Xk & B} 

n^oo fl 

which denotes the relative frequency of occurences of Xk in the set B. 


Theorem 1.4.2 Let Xk be as in equation (1.1) and assume Xq G F. Then 




for any Borel-set B C 


Idea of proof: Let O be a feasible open set for the Open Set Condition. Then O D F, 
z/(0) = 1, and for any finite sequence ii,...,ik in the set : = 

Si^o...o Si^{0) will satisfy = [1^= r® which is also the relative frequency 

of the occurences of the string ii,... ,ik in the infinite sequence generated by the R- 
So u is invariant in the sense of |Hut811 p.l6], and thus of the same form as ja in 
theorem 11.3.71 □ 


Now it is clear that the sequence xioi, a;io 2 , • • • starting at 101 and proceeding to 


infinity will enter all pixels of F\ Lemma 1.4.1 proves that (almost) the same pixels 


will be run through by regardless of whether or not Xq G F. Corollary 1.3.8 and 


theorem 1.4.2 say that F equals the support of i/, and since every pixel of F contains 
at least one point of F each one of these has a positive probability of being marked 
black while iterating Xk- 

The internet provides a great variety of software which is using this principle to 
generate self-similar and also self-affine sets. For the generation of the pictures in 
section 


6.1, the program “Fractal Explorer” |fra] has been used. 












Chapter 2 


Fractal Curvature 


Curvature Measures are well understood for sets such as differentiable manifolds or, 
more generally, for sets of positive reach |Fed59j . The hrst subsection summarizes 
the main results for curvature measures in classical geometry, but in a slightly less 
general setting that also serves our purposes: We are only going to consider sets that 
lie in the convex ring TZ of hnite unions of compact convex sets. Then we describe 
how Fractal Curvatures can be dehned for sets F whose ^-parallel sets are in TZ for 
arbitrarily small e, and hnally formulas for the fractal curvatures are given if these 
sets are self-similar. 


2.1 Curvature measures on the convex ring 


Steiner-formula. Most of the statements from this subsection are proved in |SW92j . 
We write IK for the system of all compact and convex subsets of M'^. First of all, recall 
the dehnition of the e-parallel set F^ for subsets F C with the euclidean metric 


d(-, •) (theorem 1.1.3), and assume that iF G K. Then the volume of the parallel sets 
of K is given by the Steiner-formula: 

d 

CAK.) = Y.'^'‘->Ki.,C,(K) ( 2 . 1 ) 

j=0 


Here, Cd denotes the d-dimensional Lebesgue-measure, and Kj is the volume of the 
j-dimensional unit ball. The numbers Cj{K) that are determining the coefficients of 
the polynomial in £ on the right are called intrinsic volumes. In the 3-dimensional 
case, these have the following interpretations: C^{K) corresponds to the volume of 
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K, C 2 {K) is the surface area of K, Ci{K) equals (up to a multiplicative constant) 
the mean breadth of K, and Cq{K) equals 1 (being the Euler-characteristic of the 
non-empty convex set K). 


Local Steiner- formula. For x G let pk{x) be the point in K that is nearest 
to X. There is only one such point since it' G IK. Also note that pk is (Lipschitz-) 
continuous. Write B for the system of Borel-subsets of M'^. Now, for any A & B dehne 
the local parallel set of K at A a^ 

Ue{K,A)=pfl{A)nK,. 


For hxed K, 


p,{K,-)-.= C,{Ue{K,-)) 


is a Borel-measure concentrated on K^. If both it' G IK and A E B are hxed, it can 
be shown (see |SW92j . Satz 2.3.3) that Pe{K,A) is a polynomial in e, 


Pe(K,A) 


C,(if, A), 

j=0 


( 2 . 2 ) 


and that the coefficients Cj {K, A) are hnite Borel-measures for each hxed iK G IK. 


Equation (2.2) is known as the local Steiner-formula. The term “local” arises from 


the observation that the ordinary Steiner-formula (2.1) arises as a special case of the 
local Steiner-formula (2.2), namely for A = M'^. This also shows C'j(it', M'^) = Cj{K), 
i.e. the intrinsic volumes are the total masses of the curvature measures Cj{K, •), and 
will be called j-th total curvatures later. 


Curvature measures. For each j G {0,..., d}, Cj{K, •) is called the j-th curvature 
measure of K. We note some important properties: 

1. Cj^K^Mf) = Cj{K), i.e. the intrinsic volumes are equal to the total curvatures. 

2. Cd{K, •) is supported by K, and for j = 0,..., d — 1 Cj{K, •) is supported by 
the boundary dK of K. 

^Ue{K,A) can be visualized as iG n A together with all line segments of length £ emanating 
orthogonally outwards starting from all y G A n OK. At cusps, the orthogonal directions are 
comprised in the normal cone. 
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3. Curvature measures are motion-covariant, i.e. given a rotation or translation g, 
one has Cj{gK, gA) = Cj{K,A). 

4. Ck{-, ■) is homogeneous of degree k: For A > 0, Ck{\K, XB) = X^Ck{K, B). 

5. Curvature measures are additive: Cj{K,^) = 0, and for K, M E K. such that 
K U M E K one has 


Cj{K U M, •) = Cj{K, •) + Cj{M, •) - Cj{K nM,-). (2.3) 


Note that in the last property, we were assuming KUM E M. because curvature mea¬ 
sures are initially dehned only for compact convex sets. However, via the inclusion- 
exclusion principle 

*=1 

one can formally assign a map from H to M for each hnite union IJ™ ^ Ki of sets 
Ki E K, as the intersections Hier right-hand side are members of K. It 

turns out that for different representations of the same union, [JT=1 K^ = u£i K[ , 
the maps Cj{\J^^]^Ki,-) = iF', •) coincide ( |SW921 Satz 2.4.2]). This means 

that Cj{F, ■) is well-dehned for any F E 7Z, where 71 is the ring of polyconvex sets, 
i.e. hnite unions of sets in K. 


m 

Cj(U *'..■)= E 


n^i.' (2 A) 


The maps Cj{F, •) turn out to satisfy properties 2,3,4 and 5. However, they are non¬ 
negative Borel-measures only for j = d and j = d — 1; for j E {0,..., d — 2} they are 
signed Borel-measures in general, i.e. they can have a non-zero negative variatonal 


part. For a dehnition see p,87 of the appendix. 


2.2 Rescaled Curvature Measures 

Fractal sets are in general neither polyconvex nor of positive reach, so the classical 
dehnition of curvature does not apply. However, a great amount of geometric in¬ 
formation about a fractal set F can be encoded in the behaviour of the curvatures 
Cj{Fs,-) of its parallel sets as £ 0, provided of course that Cj{F^) is dehned. In 

this section and the next, a summary of relevant results taken out of Stehen Winter’s 
dissertation |Win05j will be given. 
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To make sure that Cj{F^, ■) is defined, we assume that for every compact F there exists 
an ^0 > 0 such that F^ & TZ for all £ G (0,eo]; and we say that F has polyconvex 
parallel sets. More generality could be achieved by assuming that F^ is a finite union 
of sets of positive reach. However, this weaker assumption has not led to any deeper 
insights yet, so we stick to the setting F^ ETZto simplify the notation, knowing that 
the results can be easily extended to the positive-reach-setting. 

Consider now the expressions U'Cj{F^, A) where t G M and their behaviour as e 0. If 
t is chosen too big, they are likely to tend to 0, and for a too small t, they will behave 
eratically and probably be unbounded]^ The interesting values of t are the values at 
which e^Cj{F^, A) switches from being unbounded to being bounded. It seems wise to 
avoid cases in which e^Cj{F^, A) stays bounded but where the positive and negative 
parts e^C^{F^, A) and e^C~ {F^, A) are actually unbounded; thus the focus will be 
on e^Cj°'^{Fi;, A). To further simplify things, only the global scaling constant will be 
studied, i.e. only the case A = W^. Recall that we write Cj{F,'R'^) = Cj{F) for the 
j-th total curvature. 

Definition 2.2.1 (scaling exponents) The (upper) k-th curvature scaling exponent 
Sk and the lower k-th curvature scaling exponent Sk of a subset F C are respectively 


defined as 



Let h > 0. With this definition, one can be sure that for every A ^ B and every 
symbol • G ,nar} 


sup A) : e G (0,1]} < oo. 


and that there exists an H G H and a symbol • G {+, —} such that 


sup {e:-'' ^Cl{Fs,A): £G(0, l]}=oo. 


^Note the analogy to SNg and r® ^C'^{Ar) for the Box-counts in lemma 1.1.2 and for the 


.s — d rd 


Minkowski-content in definition 1.1.4 
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Note that an alternative way to write down definition 2.2.1 (also compare defini¬ 
tion 1.1.1 and lemma 1.1.2) is 


Sk(F) 

^{F) 


lim snp 

e4,0 

lim inf 

e4-0 


logcrjF,) 

-loge 

logCrjF,) 

-loge 




In general, Sk{F) < Sk{F), bnt for many e.g. self-similar sets eqnality prevails. 
We can now state the central definition: 


Definition 2.2.2 (fractal curvature) Let F C &e compact, and assume that 
there exists an Eq > 0 such that ^ TZ for all e G (0, eo] ■ If limits 


Cl{F) := limE'^'°Ck{F^) and 
£4.0 

(2.5) 

Cf-IF) := lim£«Cr(n) 

(2.6) 


exist, then they are called the k-th fractal (total) curvature of F and the k-th fractal 
(total) variational curvature of F. 

Unfortunately, this limit does not exist in many cases. It often happens that the total 
curvatures exhibit the growth behaviour Ck{Fs) ~ e®'', (e 0), but the expressions 
e^'‘Ck{F^) keep oscillating as £ approaches 0. In these cases, an average limit often 
exists: 


Definition 2.2.3 (average fractal curvature) Let F be as above. If the limits 


ci(F) 


lim-- 5 -- 

(54.0 log Eq — log 0 



e^’^CkiF,) 


dE 

E 


-^f,var 


f£0 


(F) := lim --^-- / 

(5t0 log Eo - log 6 Js 


(jf 

E»‘cr{E)- 


(2.7) 

( 2 . 8 ) 


exist, then they are called the k-th average fractal (total) curvature of F and k-th 
average fractal (total) variational curvature of F. 


cl{F) can be considered as a generalization of Cl{F), because if Cl{F) exists then 
so does Cf.{F), and in that case both their values coincide. The same statement 
applies to cl'^°‘^{F) and Cl’^^^{F). Note that for fc G {d — 1, d} the measures Ck{Ff) 
and Cf°''^{F^) coincide and so do their averages, so that in part II where d = 2 we 
only need to worry about Cq“^(-). 
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Consistency with classical curvature. Finally, it should be noted that the defi¬ 


nition of fractal curvature Cl {F) is an extension of the definition of the total curvature 


Cj{F) from section 2.1 If F E TZ and Ck{F) ^ 0, the scaling exponents all equal 


0, and in that case 


Ci{F) = lim£*'=Cfc(F,) = Cfc(F) 


which follows from the continuity of the measure Ck{-) and the convergence F^ ^ F 
in the Haussdorff-metric dni-,-)- Motion invariance and homogeneity also carry over 
to the fractal versions: 



Fractal curvature rueasures. As seen above, the rescaled expressions e^'°Ck{F^) 
lead to a suitable definition for the total fractal curvature of a compact set F. Simi¬ 
larly, one can consider the expressions e^Ck{F^, A) with suitable t to arrive at a local 
fractal curvature of F. The value t is best chosen independently of A, as the limits 
limgj^o will then more likely be Borel-measures in A if they exist. 

Defiuitiou 2.2.5 (weak couvergeuce of measures) Let p and pe be measures on 
the measurable space (X, S) for every e G (0,£o]- The are said to converge weakly 
against /i as e ff), if 



for every measurable function / : X —)■ M. In this case we write 


limpe = /i. 


Defiuitiou 2.2.6 (fractal curvature measure) If the limit 




exists, it is called k-th fractal curvature measure of F. Similarly, the limit 



is called the average fractal curvature measure of F if it exists. 
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2.3 Fractal curvature of self-similar sets 

If fractal curvatures exist, and what values they have if they exist, turns out to be a 
difficult question. A natural starting point for these questions are (deterministically) 
self-similar sets, since they are probably best-understood among all fractal sets. In 
his dissertation |Win05j . Steffen Winter showed the existence of all fractal curvature 
measures for a subclass of the self-similar sets and gave formulas for their calculation. 


The class Q^. Recall that the dehnition of fractal curvature requires that F^ G IZ^ 
for all £ in some interval (0,£o]- Now self-similar sets have the following fortunate 
property: 


Theorem 2.3.1 f |LW07| . Proposition 4.6) Let F be a self-similar set. Then 

3e > 0 : F, e > 0 : F, e n^. 


In other words: Either all parallel sets of a self-similar set are polyconvex or none. 
This means that in order to hnd out whether or not there is an interval (0,£o] on 
which Fg G is satished, we just need to check F^ G for any one £ > 0. 

In order to simplify the notation we introduce the following 

Definition 2.3.2 Denote the class of all compact self-similar subsets of 

which satisfy the open set condition and have polyconvex parallel sets by . 


Scaling exponents. For Q^sets, the scaling exponents always satisfy Sk < s — k. 
Even more can be said: 


Theorem 2.3.3 ( |Win05] . theorem 1.3.2) Let F G Q^. 

the expression e^~^Cl°‘'^{Fe) is bounded. 


Then for allk ^ {0,..., d} 


In all known examples of Q^sets with non-integer dimension, in fact Sk = s — k, 
but an according theorem has not been proved yet. It seems as if for Sk < s — k to 
happens it takes a “degenerate” fractal having non-empty interior]^ 


Winter gave no theorem guaranteeing a lower bound for Sk as universally as theo¬ 
rem 


2.3.3, but at least the following criterion: 


^An example for such a set is the c?-dimensional unit cube, which can be viewed as the union of 
2^ just-touching copies of size |. 
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Theorem 2.3.4 ( |Win05] . Theorem 1.3.8) Let F G and k G {0, Let 

Train be the smallest similarity ratio of the IFS of F. Suppose there exist constants 
Eo,/3 > 0 and a Borel set B C =: 0-eo such that 

Cr{Fe,B)>l3 

for each e G {rmin^o, ^o]- Then for all e < Eq 

e^-^Cr{Fe) > Pe^o-^rPn > 0 - 


We apply this theorem to all sample sets F from section |6.1| except the Cantor dust 
and the Koch curve (these are not members of Q'^), in order to determine their 
scaling exponents: Let conv(F) be the convex hull of F. Then for each F the interior 
O = conv(F)° of the convex hull is a feasible open set. Choose e.g. Eq = O.Old where d 
is the perimeter of F and B = Osq and see that for k G {0,1, 2} and £ G ^o], 

B) > B)\> c for some value c > 0. Thus we note the following 


Corollary 2.3.5 For all k G {0,1,2} and all sets F from section 6.1 except the 
Cantor dust and the Koch curve, 


^(F) = Sk{F) = s-k. 


Formulas for (total) fractal curvatures. There are explicit formulas for the 
fractal curvatures of Q'^-sets. These involve the curvature scaling functions: 

Definition 2.3.6 For a self-similar set F associated with the IFS S = {Si ,..., Sjy}, 
the k-th curvature scaling function is defined by 

N 

Ru{e) := Cu{F,) - Y, lio,n]{e)CumF\) 

i=l 

These, together with the renewal theorem (see e.g. |Fal971 Theorem 7.2]), are the 
tools which are making it possible to derive the following formulas: 


Theorem 2.3.7 (jWmOS], theorem 1.3.6) Let F G be the self-similar set cor¬ 
responding to the IFS S = (/i,..., /at) with similarity factors (ri,..., r^v). Then for 
k = 0,... ,d the following holds: 
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1. The k-th average fractal total curvature Cl{F) exists and is equal to 


X, = - 


^s—k—1 


V Jo 


Rk{e)de, 


(2.9) 


where rj = — ’"I log^i s = dimi^(F). 


2. If F is non-arithmetic, the k-th fractal total curvature Cl{F) exists and equals 
Xk- 


Note that for k = d, the above formula holds true even if F does not have polyconvex 
parallel sets. This has been shown by Gatzouras in |Gat00j . Recall that C^iF) is 
actually nothing else than A4.{F), the Minkowski-content of F (see definition 1.1. 4^, 


- j 

and with this in mind we are going to call C 2 {F) the average Minkowski-content of 
F. Gatzouras also showed that for self-similar sets satisfying the open set condition 
always X^ > 0, a result which we are going to refer to by saying that “F has positive 
average Minkowski-content.” 


Curvature measures. Fractal curvature measures, which are the limits described 
in dehnition 2.2.6 have the following simple representation for Q^sets: 


Theorem 2.3.8 ( |Win05] . Theorem 1.5.1) Let F G and assume Sk = s — k. 


Then the average fractal curvature measures from definition 2.2.6 exist and equal 


ci(F,-)^c'(F)FJ^. 

If F is non-arithmetic, the fractal curvature measures also exist and equal 


cl(F-)-cl(F)Fi^. 


Thus fractal curvature measures for Q^gets are just multiples of the self-similar 


measures we have already encountered in theorem 1.3.7 The fact that they are all 


constant multiples of each other is the basis of the geometric invariant characteristic 
defined in section 16.41 














26 


CHAPTER 2. FRACTAL CURVATURE 



Part II 

Image Analysis 
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Chapter 3 


Limitations imposed by 
Discreteness 


The main reason for starting a new part here is the change of setting: Space is no 
longer d-dimensional and continuous, but 2-dimensional and discrete. The necessity 
to deal with discrete space, however, does not only stem from the fact that we are 
analysing images: If one thinks of real and concrete counterparts of the abstract 
fractals from part one, these have to be composed of hnitely many atom^ just as 
images are composed of hnitely many pixels. In this context, if nature is to be modeled 
by fractal geometry, the discretization problem arises sooner or later. 


The Discretization Process. We assume a rectangular grid for the pixels, which 
we model by the set 

:= {{k6,l6) e : A:, / G Z} , 

where | is called resolution. Later it will turn out useful to endow Z^ with the 
euclidean metric d{-, ■) which is rescaled in such a way that the distance between 
horizontally and vertically neighbouring pixels is 1. 

A binary image I shall be represented by 

/ : Z5 —)■ {0,1} 

^This means that in nature, a 3D-Sierpihski gasket of diameter 1 meter can at most be iterated 
around 50 times before its smaller copies reach the size of atoms. Fractals of infinite recursion depth 
only exist in the theory of part one. 
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where the value 1 corresponds to a black pixel and 0 to a white pixel. We always 
assume that is hnite, i.e. that there are only hnitely many black pixels and 

thus images are bounded. 

Finally, we represent a bounded subset A C by the binary image I a which satishes 

Ia {{k6, IS)) = 1 ^ [kS, {k + 1)5) X [IS, (/ + 1)5) n A ^ 0, 

i.e. exactly those pixels are black whose corresponding areas have non-empty inter¬ 
section with A. 


“Lacunarity”. Unsurprisingly, the fractal dimension (Hausdorff-dimension) of a 
fractal set is far from determining its structure completely. Figure 3.1 (taken from 
|Man94j i shows the hrst construction stages of one-dimensional Cantor-like dusts, 
with iterated function systems as follows: In each line which is k steps away from 



Figure 3.1: A collection of Cantor dusts of dimension 
the starred line there are 2*^+^ similarities, all of which have ratio r = In the 
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upper half, the similarities place copies of the unit interval towards the ends of the 
interval, whereas in the lower half the similarities are spread uniformly. In every case, 
the similarity dimension is the value s for which 2^^^ x = 1, i.e. s = |. The 

size of gaps, determining the lacunarity of a fractal, is high on the top and low at the 
bottom. 

Now note how the results of the discretization process differ from the topmost to the 
bottommost fractal: Going upwards from x, at some point only the hrst and last pixel 
of the unit interval will be marked black, whereas going downwards far enough will 
result in every pixel of the unit interval being marked black. Image analysers will not 
be able to differenciate between the Cantor-dust and 2 singletons resp. the whole unit 
interval. The obvious effect is that high lacunarity will result in an underestimation 
and low lacunarity in an overestimation of the fractal dimension. 


“Suitable” binary images. The degree of fractality a binary image can convey 
is strongly limited by the length of chains of still visible iterates: If one can only see 
one smaller copy of the image in itself, almost all of the fractality has been lost in 
discretization, whereas not so much is lost if one can see e.g. a smaller copy of a Smaller copy 
of a smaller copy of a smaller copy. Especially estimators based on the regression on the 
logarithmic scal^ depend on this sort of information. If the smaller copies are only 
supposed to be statistically similar to the bigger copy, additionally to a high recursion 
depth a large number of images is advantageous. 


In this sense and with the above notion of lacunarity in mind, it seems that fractals 
are more suitable for image analysis the greater their similarity ratios are. It also 
seems that by excluding fractals of very low similarity ratios around one also avoids 
cases of extreme lacunarity as on the far ends of x in hgure |3.1[ The sample images 


from section 6.1 were chosen in this spirit. 


Various definitions of lacnnarity. The notion of lacunarity arose from the need 
to characterize fractals beyond their dimension. As this concept is far too complicated 
to be expressed by a mere number, there is still no uniformly accepted definition of 
lacunarity. 

In 1994 f |Man94j ). B. Mandelbrot speaks of “shell-lacunarity” which is effectively 
the same as (the reciprocal value of) Minkowski-content. D. Gatzouras took this 

^These are e.g. the box-counting algorithm, the sausage method and the method proposed in 
chapter [5| 
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as a motivation to prove in 1999 ( |GatOO] ) that self-similar sets satisfying the OSC 
always have a well-dehned average Minkowski-content (regardless of the parallel sets 
being polyconvex or not). With formulas from S. Winter’s dissertation, the author 
calculated the average Minkowski-contents C({F) for the sets F from hgure 


3.1 


The somewhat surprising result is that at * C( (F) is minimal; preceding downwards 
C({F) tends to -|-oo; preceding upwards it increases hrst and from A; = 3 on it 
decreases and approaches the value 8.65.... Thus the inverse average Minkowski- 
content does not describe the lacunarity effect from hgure ^ ideally, as one would 
expect the starred dust to have lower lacunarity than the dust above it. However, it 
is consistent in that it is correctly monotonuous when there is no transition between 
two models (as in the step from * to the line above). 


In the same work, B. Mandelbrot mentions “gap-lacunarity”, which is (the reciprocal 
value of) Cq {F), and he has already mentioned it to be a constant multiple of C( {F), 
without even the assumption of self-similarity of the dust on the real line. (The 
author does not know if this idea has been followed anywhere in the literature, in 
spaces of any dimension.) 


In physics, lacunarity is preferably determined via the “gliding box algorithm” (see 
e.g. |PGH'*~96] ): Given an object M and a box B(r) of size r, the mass s of MnB{r) is 
measured as B is gliding through space (according to a never specihed but seemingly 
uniform distribution). Q{s,r) is the then obtained probability distribution of the 
masses r according to the box size s. The hrst two moments are denoted by Zi{r) 
and Z 2 {r), respectively, and 


A(r) := 


Ur) 


(Ur))^ 

is dehned as the lacunarity for the box-size r. Plots of A(r) against r on a logarithmic 
scale can be characteristic for certain random fractals. 


For sets as e.g. percolation clusters, the correlation of mass in different directions 
depicted by angular sectors has been examined by Mandelbrot et al, and it has been 
found that “antipodal correlation” can serve as a measure of lacunarity, too. 









Chapter 4 


A Review of Dimension Estimation 
Methods 


In order to have a reference frame for the estimators s and Ffc of the method proposed 
in chapter the author searched for methods of the estimation of dimension and 
of the (average) Minkowski-content in the literature. A large amount of literature 
has been found on dimension estimates, but nothing on the estimate of Minkowski- 
content. The reason for this probably is that, in natural sciences, the latter value 
is devoid of any meaning. However, in section 6.4, measurements of the (average) 
Minkowski-content will be used as a normalization factor for the 0th and 1st fractal 
curvatures, thus providing a new set of computable geometric invariants, at least for 
Q'^-sets. 


Of the vast literature on dimension estimates the author ruled out the part which 
deals with dimension estimates of graphs of continuous functions, or more generally 
with estimates of the fractal index a of a stochastic process. It seems like drawing 
a binary image of a graph given a time series of data does not seem apropriate for 
estimates of the fractal index, and in recent literature estimators for a are all based 
on the data themselves and not on a binary image (see e.g. |Sei05j h Thus in this 
chapter only fractal dimension estimation methods will be reviewed. 
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4.1 The Box-Counting Method 


Recall definition 1.1.1 from chapter]^ If the box-counting dimension dims(F) =: s 


of a fractal F exists, Ns has the growth behaviour 


Ns ~ C6-^ 


(4.1) 


where C is some positive constant. Taking logarithms, 


\ogNs ^ iogC + slog(l/5), 


(4.2) 


SO on logarithmic paper, the Box-counting-dimension appears as the slope of an 
asymptote. All box-counting algorithms exploit this behaviour of Ns, and a typi¬ 
cal algorithm runs as follows: 

1. Set the box-size in pixels <5 G N to a maximum value 5max- 

2. Count the number Ns of disjoint boxes of size 6 that intersect the set F in 
question. 

3. Record Ns and decrease 5. If 5 go back to step 2. 

4. Finally, plot logA^"^ against — log 5, and £t a line to the data using the least 
squares method. Use the slope of the line as an estimator for s = dimB(F). 

Typically, 5max is somewhere around one fourth to one third of the image diameter, 
whereas 5min ranges from one to four pixels. In step 3, 5 is decreased in such a way that 
the points in the logarithmic plot are being close to equidistant, i.e. 5new = L^oid/aJ 
for some constant a > 1, e.g. a = 1.4. 

Note that in step the union of all disjoint boxes that intersect F will certainly 
cover F. However, this union will not necessarily be made up of the minimal number 
of boxes possible, as would be required by the dehnition of Ns- Thus algorithms 
achieve more accurate (and usually higher) estimates if, each time step is run 
through, the box-grid is shifted several times until the minimum covering number is 
found. Finding this minimal cover is incorporated in the box-counting algorithm of 
“FracLac” (E ar02j ). a plug-in of the open-source image analysis software “Imaged” 
|imaj which has been used for the box-counting estimates from chapter 







4.2. ESTIMATING LOCAL DIMENSION 


35 


4.2 Estimating Local Dimension 


As pointed out in section 1.2, Borel-sets can be generalized by measures which are 
supported by these sets, and in this setting the notion of local dimension arises quite 
naturally. In dynamical systems, an attractor can be assigned the measure of the 
relative frequency 

/i(i?) := lim —: k < n,Xk & B} 

71^00 fl 


of a particle with coordinates x being inside a Borel-set B at discrete-time measure¬ 
ments k E N. If the attractor is “strange”, typically the local dimension of this 
measure will vary throughout the set and will not be an integer. 


In |CD89] Cutler and Dawson describe statistical methods of estimating the local 
dimension of a measure on the basis of points Xk sampled according to this measure. 
An important application of their method that they had in mind probably was to 
describe the geometrical properties of an attractor. 


The estimator. Cutler and Dawson consider the nearest neighbour statistics 

Pnix) := min IIX—x|| 

l<i<n 

where x G supp/r is a point of the support of the measure fi and the Xi are independent 
random variables sampled according to the law p. They showed that if the local 
dimension is “simple at a:”, i.e. if 

dimp,(a;) = dim/i(a;) = a{x), 


theiI3 


1 ._ ^Og{bpn{x)) a.s.^ 1 

OKA) ■ 1 ^ / \ • 

—alogn a[x) 


The positive parameters a and b are not affecting the limit of ln{x), however they 
have a biasing influence on /„(x) as an estimator for a{x). Cutler and Dawson suggest 
the values a = 1 and 6 = 2 as for a uniform distribution on an interval this choice is 
neutralizing the bias term with the lowest order in n. 


^With the conventions logO = —oo, t = -|-oo, ^ = 0 and admitting oo as a valid limit, this 
statement is also valid for x outside of the support of In that case, a{x) = +oo; if x is an atom 
of /i, then a(x) = 0. 
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Binary images. If the binary image is representing an s-set F (i.e. 0 < < 

oo), then the simplest s-dimensional measure that can be assigned to F is maybe 


M-) 


W{-nF) 
n^{F) ’ 


which spreads the unit mass uniformly over F. In order to be able to utilize the 
above method, we have to find a discrete analogon /i of /i for the binary image; to 
this end, we spread the unit mass uniformly over the black pixels. Then we recall 
from section 1.4| that we assume that the black pixels are exactly those which contain 
at least one point of the support of /i. But since this point could be located anywhere 
in the pixel, we spread the mass of each black pixel evenly across its square. The 
author believes this to be the best way of discretizing the measure fi from above. 


The algorithm. Part of the work on this thesis was the creation of the java-class 
LocalDimension in GeoStoch ( |Dep07| ) that is implementing the above method for 
binary images. It calculates a histogram of estimated local dimensions of the measure 
jl at randomly picked points. 

The algorithm implemented in the above mentioned java-class runs as follows: 


1. Sample m points (the “test-points”) according to jl in the following 

way: Record the coordinates of the black pixels in one big pixel-array B. Then 
pick one pixel according to a uniform distribution on (l,length(B)) and add a 
uniformly distributed number between 0 and 1 to both coordinates. 

2. Create an array L of length m containing the nearest-neighbour statistics l{yi) 
for each i = 1,... ,m, and initialize all entries L[i] to have the value lo{yi) = d 
which is the diameter of the image. 

3. As j ranges from 1 to n, sample the point Xj according to jl as in step 1 and 
update the array L via L [i] = min{L [i], d{xj, yi)} for each 1 < i < m. 

4. Draw a histogram for the values a{yi) := ln{yi)~^, I < i < m, where ln{yi) = 

L[i]. 

5. Choose either the modal value of the histogram or the arithmetic mean of the 
a{yi) to be an estimate for the global Hausdorff-dimension of F. 


Figure 4.1 shows a typical histogram returned by the above algorithm for the self¬ 
similar triangle A from section 6.1[ 
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dimension 


Figure 4.1: A typical histogram of local dimensions of the self-similar triangle A. 
The measured local dimensions at 20000 test points has been represented by columns 
of width 0.0025. The image has 383313 black pixels, and 300000 sample points were 
used. The average of the dimension estimates lies at 1.578, close to the real value 

1.588... 


Limitations. Note that the histogram is far from its desired form, which would 
be one single column at the fractal dimension s ~ 1.588. However, since the local 
dimension of the measure /t actually equals 2 everywhere on its support, we are glad 
to achieve this result at all. The number m of test points can be increased almost 
arbitrarily, meaning that the histogram can actually approximate the distribution of 
the ln{y)~^ as y varies over the support of fi. However, the number n of sample points 
cannot be increased arbitrarily. If it increases beyond the number of black pixels, the 
pixel-scale local dimension character of fi begins to shine through, and dimension 
estimates are approaching 2. The best value for n seems to be somewhere around 
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80% of the number of black pixels. 

For self-similar pictures which have a constant local dimension, the arithmetic mean 
of the values a{yi) yielded good estimates for the fractal dimension, except for the 
Koch-curve whose dimension was overestimated by 0.09, and the Tripet, where the 
value was 0.10 too low. 


Finally, note that measures, as opposed to sets, can be represented more accurately 
if grey-values are admitted for the pixels, meaning that a pixel can be set darker the 
more mass the measure has at this pixel, whereas a set either intersects a pixel or 
not. If in the construction procedure of the self-similar images which is using the 
self-similar measure from section 1.4| the relative frequency of entries of points in a 
pixel is measured and the grey value is calculated accordingly, it might be that a 
refined discretized measure jl will yield better dimension estimates. 


4.3 The Sausage Method 


Recall the Minkowski-representation of the Box-counting-dimension, theorem 1.1.3 
If dimg(F) exists, the shrinking rate of the Lebesgue-volume of the parallel sets 
C^{Fr) for r —)■ 0 is directly connected to the Box-counting-dimension dimB(F) via 


where AA.{F) is the Minkowski-content of F. As before (see equations 4.1 and 4.2), 
on a logarithmic scale the growth behaviour is linear. 


\ogC^{Fr) ~ log A1(F) + {d- dimB(F)) logr. 


which is the basis of the following 


Algorithm. Now proceed according to the following steps: 

1. Calculate the distance-transform of the binary image I. This is an image 
that records the distances to the nearest black pixels: 

D/ : P ^ M 

p ^ d{p,r^ {{!})) 
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Figure 4.2: A plot of the areas of the dilated binary image against the dilation radii. 


The underlying fractal set is the Sierpihski Tree from section 6.1.1 


2. For radii r that are equidistantly distributed on the logarithmic scale (the max¬ 
imal and minimal radius can be chosen the same way as S^ax and dmin in the 
Box-counting method) calculate the areas of the r-dilated binary images 

Ir:=Dj^([0,r]). 


3. Plot these areas against r on logarithmic paper and £t a least squares line to the 
data points. The slope serves as an estimate for dimB(F) — 2 (see figure 4.2). 


This procedure has been implemented by the author in the class Curvature2D of the 
GeoStoch-library as the special case 

(useEuler, useBdlength, useArea) = (false, false, true), 


for a description see section 5.2 
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It should be noted how in hgure |4^ stairs appear for small radii. In fact, since the 
dilated images R are all discrete, their area is piecewise constant as r decreases. On 
a logarithmic scale this effect is more prominent for small radii than for big radii, as 
figure |4.3| illustrates. This means that, whenever sparsely distributed dilation radii 
are used (e.g. to reduce the run-time of the algorithm), care should be taken as not 
to use those radii which are close to a jump of the area. A good strategy is to use 
the optimal-area radii, i.e. those radii for which the discretized disks Bdiscrete{o,r) 
centered at the pixel in the origin o have an area which is exactly equal to vrr^. With 
the discretization mentioned above, i.e. 

Bdiscrete{o,r) := {p E P : \p - o\ < r} , 

where p and o are the coordinates of the center of two pixels, these radii can be found 


as the horizontal intersections of the two graphs in hgure |4.3[ whereas the jumps 
occur at the vertical intersections. 


The hrst optimal-area radii (rounded to 4 digits) are 

0.5642, 1.262, 1.696, 2.585, 3.432, 3.785, 4.406, 4.687, 5.322,... 
the exact values of which are 



corresponding to the areas 

1,5,9,21,37,45,61,69,89,... 

More optimal area radii of size up to 400 are listed in the class OptimalRadii of 
the GeoStoch library ( |Dep07| ). These numbers have been found numerically with 
Matlab as the zeros of the function 

V{r) = C\Bd^screte{o,T)) - C\B{o,t)). 

Also see the paragaph on the sausage method by Stoyan and Stoyan |SS94j . who 
suggest admitting corner points of pixels as possible centers of discrete circles. 


4.4 Further Methods 

Variations of Box-counting. For images of deterministic fractals, variations of 
the box-counting method have been proposed: Sandau and Kurz suggest the ’’Ex¬ 
tended Counting Method” f |SK97j 1. which is depicting the highest dimension mea¬ 
surable by box-counting if the position of the coarsest grid is varied. Martinez-Lopez, 
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Cabrerizo-Vilchez and Hidalgo-Alvarez f |MLCVHA01] ) suggest replacing the number 
Ns{E) of boxes of size 6 by the number 

OO 

j=i 

where are the (just touching) boxes of size 6 used for covering F, | ■ | is the 
diameter of a set and s is the (unknown) Hausdorff-dimension. An initial value Sq is 
chosen, and iteratively Sj+i is computed as the least squares fit of a line to the plot 
of Ng^,{F) versus — log5 on a logarithmic scale. It was shown empirically that s* 
converges very quickly. 

Stochastic processes of fractal index. A large class of random fractals for which 
dimension estimates are important is generated by stochastic processes. As an ex¬ 
ample, rough surfaces are often modeled by graphs of stationary gaussian processes: 
These are families of real-valued random variables satisfying 

At ~ Ar(0, a^) Vf e R^. 

Adler f |Adl81j l showed that if the covariance function 

7(t) := cov{Xo,Xt) 


Af> n F 
5 


satisfies 

j(t) = j{o)-c\tr+o{\tn (f^o) 

for some 0 < a < 2 and c > 0, then the graph 

Gx := {{t,Xt) e R^+^ -.teR^} 

of At has Hausdorff-dimension d — f almost surely. Hall and Roy ( |HR94j ) extended 
this result to processes of the form g{Xt) for certain smooth g and At is as above. 

Other important examples are the paths of fractional Brownian motion (fBm) with 
Hurst-parameter H G (0,1): This is a (non-stationary) centered gaussian process Bt, 

i.e. 

Bt ~ A/'(0,(Tt) 

with covariance function 

7(s, t) = cov{Bt, Bs) = E{BtBs) = ^ - \t - s\^^) . 
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It is not hard to see that the two processes {a~^Bat, t > 0} and {Bt, t > 0} have the 
same distribution, and thus Bt is said to be statistically self-similar. The paths of Bt 
are fractal sets of Hausdorff-dimension 1 + H. 

Thus the knowledge of the dimension of these random fractals is equivalent to the 
’’fractal index” a resp. H of the according processes, and in the last 20 years a big 
effort has been put into constructing estimators of the fractal index. A good summary 
can be found in the recent PhD-thesis of Sei Tomonari ( |Sein5] ). 
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-log(radius) 


Figure 4.3: The areas of the discrete and continuous disk, as a function of the radius. 
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Chapter 5 


Image Analysis with Fractal 
Curvature 


In the following we are going to propose a dimension estimation method that exploits 


the connection between the scaling exponents Sk from dehnition 2.2.1 and the dimen¬ 
sion s := dim(F) of a fractal set F. Similarly to the sausage method, we are going 
to consider the parallel sets F^, but additionally to the area C 2 {F^) we are going to 
measure the other (total) curvatures Ci{F^) and C'o(Fj as well (i.e. half the bound¬ 
ary length and the euler number), and perform a linear regression simultaneously on 
all three data sets. It appears that the quality of the measurements of the Ck{F^) 
is decreasing with k, but in several cases, especially for non-arithmetic F, the addi¬ 
tional data increase the quality and robustness of the dimension estimate. We assume 
F E throughout this chapter, and that the scaling exponents satisfy Sk = s — k. 
As noted in section 2^, this seems to be the case whenever F E and dimF ^ Z, 


and it has been verihed for the sample sets of section 6.1 in corollary 2.3.5 


5.1 The Regression Model 

The logarithmic scale. Since iterated function systems shrink a set E with the 
ratio Ti every time a similarity Si is applied, the size of the sets Sk o ... o Si{E) 
decays exponentially with k. Therefore logarithmic rescaling is useful to describe 
growth phenomena for these fractals. Now let us express the growth behaviour of the 
curvatures of the parallel sets on a logarithmic scale. First, recall that two functions 
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f and g : D are said to be asymptotic to each other for x —)■ a G -D iQ 

fix) 


gix) 


= 1 + o(l), X —)■ a, 


provided of course g{x) 7 ^ 0 in some open set containing a. 
Now the renewal theorem tells us that 


(5.1) 


where pk : (0,1] —)■ M is either constant (if F is non-arithmetic) or periodic of mul¬ 
tiplicative period e~^ (if F is h-arithmetic), i.e. pkie~^e) = Pki^)- Note that the 
expression on the left of (5.1) is positive for all £ > 0 and all fc, so is well-dehned. 


For brevity we introduce the notation 

1 


/ := lim 


/(^)— = lim [ fie ^)dt 


<54.0 — log 6 Jg ' E ^ ± j Q 


( 6 . 2 ) 


T-s-oo T 

if the limits exist for a function / : (0,1] —)■ M, and call it the average of f. Note 
that functions that are asymptotic to each other have the same averages, so with 
gkie) := e^-^Cl^^{Ff) we have 

y^f,var, , 

Pk = 9k = (F). 

For a concise notation we decompose pk into two parts: 


Pkie) = ci’''''\F)pl{e), 


so that p^ = 1. Now (5.1) can be rewritten as 


^ c’r\F)pi(E)E- 


(e 4 . 0 ), 


(6.3) 


(5.4) 


and after taking logarithms and setting e := e ^ we hnally arrive at 


log 


/^var ( TP 
c>—xk 




logC'(,’““" (F) + logp°(e -1- sx, (x -)■ CX)). (5.5) 


Recall that Cp'^{Ff) is just the area of F^, and C'Pf^{Ff) is just half the length of 
the boundary of Fg. In general we have no way of measuring Cp^(-) with an image 
analyser, because it is not additive on IK and thus the inclusion-exclusion principle 
does not work. Sometimes, however, Cp^{F^) can be precisely calculated as it is just 
equal to 2iV(Fe) — Fo(F£), where N is the number of connected components (see p.83 
of the appendix). In case Cp^{F^) cannot be calculated, it will be excluded in the 
following regression. 

^By 0(1) we just mean a sequence converging to 0 . 
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Simultaneous regression. Let X = {xi,..., be a set of x-values. Having 
equation (5.5) in mind, introduce the following variables: 


Ukj 


(5.6) 

Dk 

:= 

(5.7) 

Qkj 

:= log plixj). 

(5.8) 


Equation (5.5) suggests that, given that the approximation 


” is good enough, for 

hxed k the Ukj all lie on the graph of a function which is the sum of an affine function 
Dk + sx and a periodic function qk{x) (of additive period h if F is h-arithmetic). 
Figure [ffiT] shows a plot for the Sierpihski Carpet of size 3000 x 3000 pixels, where the 
corresponding |/fc-values have been calculated with an image analyser. This suggests 
the following regression model: 


(6.9) 


where the Skj are random variables which are modelling discretization and measure¬ 
ment errors. Since we do not know any better, we assume (rather optimistically) that 
for each k and j they are normally distributed 5kj ~ A/'(0, cr^) for some cr^, and that 
they are mutually independent. 

The problem is now to estimate the unknown values s anc|^ Fq, ..., Dd, given some 
X-values Xj, j G {1 ,... ,m} and the corresponding d + 1-dimensional point cloud 


(yoj\ 


(Do\ 


( 


(s\ 


AoA 


= 


+ 

(lij 

+ 

s 

Xj + 

6ij 



\DdJ 


WjJ 






3 ^:= 


oi 


) • • •) ^dj) ■ Tkj Yk{xj), j 1, •••, 


m 


Note that the the points of 3^ concentrate around a one-dimensional line segment in 
whose length is related to the fractal dimension s and whose base point is given 
by (Dq, ..., Dd) (assuming xq = 0, see hgure |5.2[ ). One way of obtaining estimators 
and is to ignore the periodic components of t/kj and to ht a line to y so 
that the sum of squared distances is minimal. 

^The variables Dk turn out to be averages of the logarithms of the fc-th total variational curvatures 
of the parallel sets. By putting Cl{F) := exp(Z 3 fc) one could define another average of the fractal 
curvatures; the author did not follow this idea as he had no desire to reformulate all the results of 
Winter’s dissertation. 
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Figure 5.1: A plot of pk against x for the 2D Sierpihski Carpet. Notice that for too 
big dilation radii (i.e. r > exp(5) x < —5) CQ°''^{Fr) turns to 1, and the logarithms 
become 0. For too small dilation radii (i.e. r < 1 x > 0) the measured curvature 
Ck{Fr) is piecewise constant as the maximum pixel-resolution is reached, yielding the 
slopes 0,1,2 for j/q, pi, P 2 respectively. 

Proposition 5.1.1 For a self-similar set F G consider the simplified regression 
model 




(Do\ 


(s\ 


AoA 

yp 

= 


+ 

s 

Xj + 

6ij 



\DdJ 


w 




(5.10) 
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Figure 5.2: The point cloud y 


with variables as defined in (5.9). Then the least squares estimate {s^^\ 
for (s, Do,..., Dd) is given by 


D 


(rr^)^ 


where x ■.= — Xj, 

m ^-^0 = 1 J ^ 


ym) _ _ V _ L 

m{d + T){x‘^ — x'^) 

Df) = Vo-xs 
bf'^ = Yd-xs 

Y-.= ^Y.7^,x^,<^ndY:=iY.7^,Y- 


Calculation: Differentiation of the sum of squared errors 

d m 




m{d + 1) 




i=0 j=l 
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yields the critical point ..., which has to satisfy 


(<i+i)iyM + *£.<”» + ... + = iEto El’ll 




+ ^ 


(m) 


Vn 


xs 


(m) 


+ 


D 


(m) 


Yd, 


a linear system of equations, the solution of which yields the above results. 


□ 


If the periodic components pk of (5.9) are ignored, one can nevertheless arrive at a 
reasonable estimate s for s: 


Proposition 5.1.2 Assume that the regression model (5.9) holds true and let X = 


{xi,X 2 ,...} be such that Xj —)■ oo, and write for the least-square estimator for s 
corresponding to the first m data points. Then, with probability 1, 


;(m) 


—>■ s, (m —)■ cx)). 


Proof: This is a consequence of 


Yu 


kj 


—)■ S 


Xi 


almost surely for all fc = 0,1,2 and j ^ oo. 


□ 


However, for a choice of X as above, the least-squares estimates that correspond 
to the first m data points converge to Dk only for non-arithmetic sets F in general: 


Proposition 5.1.3 Let F E Q'^ be a non-arithmetic self similar set, and assume the 


regression model (5.9) holds true. Let and be the least squares estimators 


for the simplified regression model (5.10) as defined in Prop 5.1.1, associated to the 
first m data points. Then 




{m -E oo) for all k E {D,..., d} with probability 1. 


Proof: Since F is non-arithmetic, the normalized periodic functions p^k{s) from (5.3) 


are constant and equal 1, and thus the simplihed regression model (5.10) is equivalent 


to the regression model (5.9). Using standard results on linear regression one sees 
that the estimators converge almost surely. □ 
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Proposition 5.1.4 Let F & be an h-arithmetic self similar set, and let and 
be as above. Furthermore, assume that the set of x-values X = {xi,X 2 , ■ ■ ■} is 
uniformly distributed mod h and that it is converging to +cx). Then 

ess sup > ess inf 

with respect to the Lebesgue-measure on M implies 

limsup-D^™^ < Dk for a// A; G { 0 ,, d} with probability 1, 

m^oo 

i.e. is asymptotically biased for each k. 


A typical set of x-values with the above property is A = {aj : a > 0, j G N, | ^ Q}. 

Idea of proof: If F is h-arithmetic, the estimators and f)^^\ k & {0,... ,d} are 
minimizers of the sum of squared deviations 




(m)x 


inf 


q,E)eiRxR'«+i m{d 1) 


d m 


i=0 j=l 


Now, as m —)■ oo, we can change the order of the infimum and the limit, and 




inf 


^ a ^ m 


{t,E)mxR<i+i d-\-1 '^m- 5.00 m 


.. d .. m 

= inf -—- lim — 'S^{Yij — {Ei-L sxOy 

i=0 j=l 

^ d ^ m 

= inf lim —y^(Di - Ei-h Qij + S, 

EgR‘‘+i d +1 m^co m 

i=0 i=l 


1]) 


= inf 


d + 1 
1 


((A - Eff + + qf + 2(A - A)g. 


i=0 

d 


= inf 


d + 1 


i=0 


(((A - A) + QiY + ~ + 


a 


a unique inhmum which is attained if Ei = Di + qi for all i G {0,... ,d}. Note 
that in the second equality without loss of generality we assumed t = s, a.s for 
t y s the limit is 00 . In the fourth equality we have put g* := logp^{e~^)dx and 



52 


CHAPTER 5. IMAGE ANALYSIS WITH ERACTAL CURVATURE 


/o^ I log Pi* ^)\^dx and used the ergodic theorem. Thus converges to 

Di + Qi. With the Jensen-inequality, one sees that 

j-h 

\ogp^^{e~^)dx < log / p^{e~^)dx = 0 

Jo 

where the inequality is strict because the logarithm is a strictly concave function and 
because we have assumed not to be constant. □ 



Estimating Fractal Curvature. As we have seen, for arithmetic sets F & 
the estimator will systematically underestimate logCf(F), even if we assume 
that X can tend to cxd , i.e. that we can represent e-parallel sets for arbitrarily small 
e. This is essentially due to oscillations of the rescaled curvatures that average out 
negatively under the concave function log. But even for non-arithmetic sets, plots 
of rescaled curvatures may also show oscillations, since these might only vanish for 
inhnitely large x. Therefore, a reasonable estimator for the fractal curvature Cj^{F) 
must smoothe out these oscillations, and this can be achieved via averaging as in 
: Recall that Ci.{F) equals the average of the rescaled curvature function, i.e. 

1 

cl{F) = lim — / exp{—x{s — k))Ck{Fe-^)dx. (5-11) 

T^oo T Jq 

In practice, however, we can only estimate this integral, as Qkis) can only be deter¬ 
mined for hnitely many values of e, and because we can only use an estimate s for 
s in the above formula. The integral will be estimated by the area of a histogram 
whose j-th column is centered at Xj and has height equal to Hkj for j G {1,..., m}. 
To this end we dehne the endpoints of the histogram to be 



fg •— 

^ __ Xj+ Xj+i 

h - 2 

tm • T ® 


, 1 < j <m - 1, 


where a is half a typical stepwidth between the Xj. Thus the limit in (5.11) can be 
estimated by0 

1 


r(™) ._ 
k ■— 


tm Iq 


^exp{-Xj{s - k))Ck{F^-^j){tj - (5.12) 


i=i 


^Equation (5.121 seems to be computationally unstable; instead, the author used the form 

{CkiF-.Ry 


1 


to 


rn / 

^ sgn(C'fc(E-X,)) exp f -sx^ -f log 
i=i ^ 


(tj b-i) 
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As seen before, \Cl{F)\ is underestimated by , under the assumption of course 
that the regression model above is valid. However, note that 

< Tk 


whenever both estimators have the same underlying data set. The size of their relative 
difference is inversely related to the quality of the £t of a line to the according data 
set: Usually the fit of a line to j/ 2 . is better than to yo., meaning that for k = 2 the 
relative difference between the above two estimators lies somewhere between 1% and 
, whereas for A; = 0 it typically ranges from 5% to 12%. 

is slightly more accurate than 


Applied to the sample images from section 


6.1 




5.2 Algorithms 

From now on we assume d = 2 until the end of this thesis, as we describe the 
implementation of algorithms for the calculation of and for /c G {0,..., d} 
which we have yet only done for 2D-images. A generalization to higher dimensions is 
straightforward, it should be noted however that for k < d — 2 the total variational 
curvatures may be hard to obtain. 

An algorithm for the estimation of s = dimF and Cj^{F) for k G {0,1,2} has been 
implemented by the author and included in the GeoStoch |Dep07| library. It runs as 
follows: 

1. Set the parameters: 

(a) Set the boolean variables useEuler, useBdlength and useArea to true 
if the corresponding data yk- for k = 0,1, 2 shall be taken into account in 
the regression. By default, all three variables are set to the values false, 
true, true respectiveljj^ 

(b) If quickEvaluate is set to true, the array of dilation radii R is chosen 
as explained in the paragraph below. If not, the radii of R will be chosen 
uniformly on the logarithmic scale, according to the parameters r_min, 
step and r_max. 

^Note that a choice of false, false, true will return a dimension estimate by the sausage 
method from above and an estimate of the average Minkowski-content. 
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2 . Calculate the distance-transform Dj of the binary image / which represents the 
fractal set F. This is an image that records the distances of each pixel to the 
nearest black pixel: 

Dj :P M 

p ^ d{p,r^ {{!})) 


3. As Tj runs through R in an increasing order, represent the parallel set by 


F,, :=D,-'(|0,r]), 


and for each k G {0,1, 2} estimate C^‘^''{Fr.). For k G {1, 2} this is done by cal¬ 
culating the k-th Minkowski-functional C'fc(Fr^)0the 0-th variational curvature 
C^°-^(^Frj) is estimated by -|- Q = 2N — Co{Fr.) where N is the number of 
connected components of F^- and Q is the number of holes of F^... (This goes 
wrong quite frequently, see p,83 of the appendix.) 

For the calculations of the Minkowski-functionals the algorithm proposed by 
Klenk, Spodarev and Schmidt (see |KSS06j ) and Guderlei f jCKM+OTj ) is em¬ 
ployed, which for each binary image calculates all of the above three functionals 
simultaneously. Now the values ykj are calculated as in equation (5.6), and Xj 
is set to — logTj. If the boolean variable brk is set to true, then the algorithm 
will jump to the next step as soon as -|- Q < 2. If not, a warning will be 
displayed, saying that the dilation radius might already have grown too big for 
Fr- to reveal any fractal structure. 

4. The estimator s from the simplihed regression model (Proposition |5.1.1 ) is 
calculated, based on each data set Hk. for which the corresponding boolean 
variable (useArea, useBdlength, useEuler) is set. 

5. The estimator F^™^ of the fc-th total fractal curvature of F is calculated via 


averaging as in (5.12). 


Note that the data set {y^j : j G {1,..., m}} should be excluded by setting useEuler 
= false if CQ°‘'^{Fr.) cannot be calculated properly, as in that case the yoj do not 
show growth behaviour suitable for regression. The break condition 

‘W + Q < 2” 

in step 1^ helps to avoid the case yoj = logO. 


®and dividing by 2 if fc = 1. 
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The choice of dilation radii. The proper choice of dilation radii is essential to 
the accuracy of the estimates s and k E {0,1, 2}. Some thought has to be given 
to the range and to the distribution of dilation radii: 


• If data for r G (0,1) are included in the regression, the author has observed 
that the slope of y 2 will be slightly underestimated, while the slopes of yi and 
j/o will be overestimated, meaning that the overall accuracy decreases. Thus 
only radii > 1 should be used. The default value for is 1.261. 

• If the radii are increased too far, TV will completely lose its fractal character. 
A reasonable upper bound for the radii can be e.g. 

max{0.06 * min{t(;, h}, 20} 


where w is the width and h is the height of the image. This is also the default 
value. 


Since the parallel sets Fr are represented by binary images Fr, the curvatures 
Ck{Fr) will only assume discrete values, and Ck{Fr) will be piecewise constant]^ 
Thus the main amount of information, encoded in yk on an interval where 
Ck{Fg-x) is constant, is extracted by the hrst x that falls into this interval; 
Further values x that fall into this interval will only have an averaging effect. It 
seems reasonable to pick radii somewhere out of the middle of these intervals, 
and not close to their ends which are the points of discontinuity of the Ck{Fr)- 


These points are given e.g. by the optimal area radii from section 4.3 


With this in mind, the author suggests two different choices for the set of dilation 
radii R, depending on whether or not the run-time of the algorithm is of importance: 


Lower run-time: The hand-picked radii quickRadii from the class OptimalRadii 
are used. These are radii yielding the optimal discretized area, starting from 
'T'l = fmin = 0.5641 and increasing in multiplicative steps of around 1.5. 

Higher precision: Choose 

R = {rj : Tj = rmmStepJ, j G N,rj < 

max } 

®Note that on a logarithmic scale, this effect is more prominent for small radii than for big radii, 
see figure [T3| Also observe this effect in figure [O] where the yk are piecewise linear as a function 
of a: = — log r, the slope being equal to k. 
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where the multiplicative stepwidth step > 1 can be chosen almost arbitrarily 
close to 1. Default values are = 1-2616, r^ax = max{0.06 * min{tn, h}, 20} 
and step = 1.05. 



Figure 5.3: A plot of i/o for the Sierpihski Carpet SC and for two different data sets. 
Both data sets can be thought of as being comprised of 4 full periods. However, 
extrapolating to inhnity one sees that the slope of £t 1 is preferable to the slope of 
£t 2. Indeed, the slopes are 1.87 (£t 1) and 2.04 (£t 2), and the dimension s of the 
Sierpihski Carpet is 1.89. With the breaking condition ‘W + Q < 2” one achieves that 
the lowest stair is not aborted, i.e. too high slopes as in £t 2 above will be avoided. 


Problems. For arithmetic self-similar sets the following problem arises: As the 
range of radii, for which the values y^j are calculated, is bounded, there will most 
probably be an incomplete period in the data set. But even if there is no incomplete 
period, the slope can severely deviate from its correct value, depending on where the 


beginning of a period is put (see hgure 5.3). The consolating thought here is that 
examples of arithmetic self-similar sets in nature arise even less often than examples 
of non-arithmetic self-similar sets. 






Chapter 6 


Results 


6.1 The Sample Images 






k 

eA * . 





DEEH 





Figure 6.1: The Sample Images 


2 

As most canonical examples of self-similar sets as e.g. the Sierpihski gasket, Sierpihski 
carpet, Cantor dust and Sierpihski tree are arithmetic in the sense of section [2731 some 
effort has been put into constructing self-similar sets that are non-arithmetic (and 
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have polyconvex parallel sets). These sets are the non-arithmetic self-similar triangle 
and square below. It turns out that, as e 0, the growth of especially is 

more uniform for non-arithmetic sets than for arithmetic sets, which allows a better 
£t in the regression analysis of chapter 


The sets from this chapter all have polyconvex parallel sets except the Cantor dust 
and the Koch curve. To see that, imagine a dilation of the set in question by a large 
e; it will be equal to the ^-parallel set of the set which is the union of the filled-in 
smaller copies, which is of course polyconvex. Applying theorem 2.3.1 one now sees 
that all parallel sets are polyconvex. 


6.1.1 Canonical Examples 

These are the Sierpihski Gasket, Sierpihski Carpet, Sierpihski Tree, Koch curve and 
the Cantor Dust. 



Figure 6.2: The Sierpihski Gasket 

Sierpihski Gasket SG. Here the Iterated Function System consists of three sim¬ 
ilarities which map (equally oriented) half-sized copies of the gasket towards its ver¬ 
tices. All three similarities have the ratio and thus the gasket is log 2-arithmetic, 
and its dimension s is the solution of 3 = 1 s = ~ 1.585.... It has poly¬ 

convex parallel sets. The 0-th variational measure Cq'^^(SGs) of the parallel sets can 
be calculated accurately by an image analyser, being equal to 1 — Co(SGs) 
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The author recalculated the fractal curvatures of SG in his dissertation f |Win05j h If 
the base has length 1, the correct values are: 

C^(SG) ^ -0.042345... 

C((SG) ^ 0.37615... 

Ci(SG) ^ 1.81 



Figure 6.3: The Sierpihski Carpet 

Sierpihski Carpet SC. The Sierpihski Carpet SC is the attractor of eight simi¬ 
larities of ratio | and thus has dimension s = ~ 1.893.... It is log 3-arithmetic 

and its parallel sets are polyconvex. 

On a scale where the base line has length 1, the average total fractal curvatures are 

Ci^(SC) ^ -0.0162 
C((SC) ^ 0.0725 
Ci(SC) ^ 1.352 

as computed by Winter, and C((SC) was rechecked by the author. 

Sierpihski Tree ST. ST is the attractor of three similarities, all of ratio A 

1 . rotation by — | and placement in upper left corner 

2 . placement in lower left corner 
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Figure 6.4: The Sierpinski Tree 

3. rotation by +| and placement in lower right corner. 

Similarly to SG, ST has dimension s = ~ 1.585... and is log 2-arithmetic. It has 

polyconvex parallel sets. But unlike the Sierpinski gasket, here we have Xo(ST) = 0 
and in fact the 0-th fractal curvature Cq(ST) = 0. Furthermore, Cq°'^(ST^} cannot 
be measured correctly, and thus in the regression process of chapter [^useEuler has 
to be set to false. 


Figure 6.5: The middle third Cantor dust. 


Cantor Dust CD. Recall the Cantor Dust in one dimension: In each step of the 
construction, the middle third of the unit intervall is erased. In the two-dimensional 
case, in each step a cross of thickness one third is erased. 
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The corresponding Iterated Function System consists of 4 similarities of ratio 
mapping CD in the upper left, upper right, lower left and lower right corner. As 
all similarities have ratio |, CD is log 3-arithmetic, and dim(CD) = cs 1.2619. 
This is an example of a self-similar set with no polyconvex parallel sets, so strictly 
speaking the theory of chapter does not apply. However, there does seem to be 
the correct scaling behaviour S 2 = s — 2, si = s — 1, and dimension estimates based 
on the regression method of chapter]^ with (useEuler, useBdlength, useArea) = 
(false, true, true) are quite exact. 



Figure 6.6: The Koch-curve 

Koch curve KC. The Koch-curve KC shares the dimension-, arithmeticity- and 
polyconvexity- properties of the Cantor dust CD, and also seems to have the same 
scaling behaviour for S 2 and si with equally good dimension estimates using the above 
parameters. It is noteworthy that the local dimension estimates are too high for this 
set. 

6.1.2 Further Examples 

Modified Sierpihski Carpet MC. The iterated function system of the modihed 
carpet consists of 8 similarities of ratio |, just as for the original carpet, so we have 
the same dimension, s = ~ 1.893..., and log3 arithmeticity. There are some 

differences in the positioning of the smaller copies (there is a copy placed in the 
middle now) and some of the smaller copies are rotated by multiples of The 
parallel sets are polyconvex. Fractal curvatures are taken from |Winn5j . 


C'o^(MC) ^ -0.014... 
C/(MC) ^ 0.0720... 
Ci (MC) 




1.344... 
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Figure 6.7; The modified Sierpinski Carpet 



(a) The Tripet (half triangle, half carpet) (b) The similarities generating the 

Tripet. Thanks go to Steffen Winter for 
this idea and picture. 

Figure 6.8: The Tripet. 

The Tripet. Like the Sierpinski Carpet and the modified Sierpiiiski Carpet, the 
tripet is the attractor resulting from 8 similarities of ratio | and thus has dimension 
s = ~ 1.893.... It is also log 3-arithmetic. The specialty abont this set is that 

standard methods seem to underestimate the dimension. 

Non-arithmetic self-similar triangle A. Consider the self-similar triangle A 
from the above figure with a right angle at the bottom left corner. Two smaller, 
positively oriented copies (green and bine) of the triangle are placed in the bottom 
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Figure 6.9: The non-arithmetic self-similar triangle A 

corners, while the similarity that maps the triangle to the upper vertice (red) is a 
scaled reflection and thus negatively oriented. A few calculations show that, in order 


for the vertices of the smaller copies of the triangle to meet as they do in hgure 6.9 
necessarily 


ac 


X = 


IP -L 


y = 


bP 


Pc 


z = 


^ 62 + A’ 62 + c2’ 

and thus the above triangle is unique up to similarity transformations. To fix the 
size of A we choose a = 0.6, 6=1 and c = 0.8, then the similarity ratios r 2 for the 
reflection and ri, for the remaining two similarities become 


25 


20 


ri = 


r2 = 


rs = 


16 

41‘ 


41’ 41’ 

The similarity dimension s which is the solution to 

r* + r2 -F rg = 1 

was computed numerically with Matlab, yielding the result 

s 1.5882... 

The set A is “properly non-arithmetic”, in the sense that 

log Ti 




log Tj 
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To see that e.g. ^ Q, suppose that 

P r ^ ^ ^ ^ 

--= - tor some p G A, q G N, 

log r 2 q 

then (using the uniqueness of prime number decompositions) 

rl = rl^ AU-P X 5“'^ x = 1 p = g = 0, 
a contradiction to g G N. 



Figure 6.10: The parallel sets of the (smaller copies of) A. 


As A is a (non-arithmetic) self-similar set with polyconvex parallel sets, its fractal 
curvatures A/(A) are well-defined. The tedious task of computing them via theo- 
2.3.7| is described in the appendix on page ^ here only the results shall be 


rem 


given. 


Ao^(A) ^ -0.023459108... 
Ci{/\) ^ 0.239312913... 
Cl{/\) ^ 1.162171558...., 
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on a scale where the hypotenuse has length 1. 



Figure 6.11: The sheared Sierpihski Gasket 

The sheared Sierpinski Gasket SSG The Sierpihski Gasket can be “sheared” 
in such a way that it has the same convex hull as the self-similar triangle from above. 
The dimension stays the same: dim(SSG) = dim(SG) = ~ 1.558.... For a 

comparison of measured total curvatures see pages and [76| 



Figure 6.12: The sets Ml, M2 and M3 
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The sets Ml, M2 and M3. The three polyconvex self-similar sets from hgure |6TL2 
are the attractors of 24 similarities of ratio I and thus all have the dimension 

D 


S = 


log 24 
log 6 


1.7737... 


Since the similarity ratios are relatively low, there are not so many iterations possible 
before the maximum resolution is reached, and naturally this means that discretiza¬ 
tion errors will be higher compared to the other sets (except the window and gate 
below which similarity ratios |). The “gaps-sizes” are decreasing from Ml to M3, 
which is reflected in increasing estimates of dimension and of curvatures. 




(a) The Window (b) The Gate 

Figure 6.13: Two versions of the Sierpihski Carpet 

The Window and Gate. The dimension of these two sets is ~ 1.8957... 

log 7 

as there are 40 similarities of ratio The author calculated their curvatures in 
calculations which were even longer than those for the set A and which are not given 
here for the sake of brevity. The results are: 

C^{W) ^ -0.0146171712902 C^{G) ^ -0.0163916537451 

C((W) ^ 0.0652764265706 C((G} ^ 0.0732007965716 

Gl{W) ^ 1.251813666054 C'|(G) ^ 1.403780236274 

On the pixelscale, i.e. where the base has length 2920, the fractal curvatures rescale 

log 40 

with 2920 , yielding 

G^{W) ^ -54228 G^iG) ^ -60811 

G( (W) ^ 242166 G( (G) ^ 271565 

G^{W) ^ 4644054 gI{G) ^ 5207828 
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It is noteworthy that the specihc fractal curvatures (i.e. 0th and 1st curvatures if the 
2nd curvatures are normalized to 1, see section 6.4) are identical up to 14 signihcant 
decimal hgures: 


So(hh) 

Si(hh) 


So(G') ^ -0.011676794 
Si(G') ^ 0.052145481. 
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39 
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Figure 6.14: Left: a non-arithmetic self-similar square. Right: its subdivision by 
similarities. The big square has base-length 100, the numbers show the base-lengths 
of the small squares. 


Non-arithmetic self-similar square □. Figure 6.14 shows a self-similar square 
that consists of eight small copies of itself like the Siepihski Carpet, but the sizes 
of the copies vary according to the numbers on the right. Computing the similarity 
dimension numerically yields the value s = 1.7937.... Note that □ is non-arithmetic, 
and that C'Q“'’(n) cannot be measured properly. 


6.2 Estimates of Dimension 

The estimator s of the fractal dimension which is based on the multiple regression 
algorithm from chapter will be compared to the already known estimators, a short 
description of which can be found in chapter The sample sets from the previous 
section will be binary images in the versions “large” (3000 x 3000 pixels) and “small” 
(250 X 250 pixels). For the small versions, cutouts of the large images have been 
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chosen, and care has been taken that the cutouts do not disappear completely in a 


hole of the big fractal. Tables 6.2 and 6.2 show comparative results. 


For the box-counting estimates we chose the program “FracLac” ( |Karn2] ). a plug¬ 
in of the open-source image analysis software “Imaged” |imaj . The local dimension 
estimates were obtained from an algorithm that the author has implemented himself 


(see section 4.2). It has been sorted into the “GeoStoch”-library ( |Dep07| ) as the 
class LocalDimension. The default parameters for this method were used in the 
below measurements: 1050 test points in which the local dimension is estimated, and 
0.8 X the number of black pixels has been used as the number of sample points for the 
nearest neighbour characteristics. Unlike in |SS94j , we do not choose the modal value 
of the histogram of dimension estimates as the estimate for the global dimension, but 
the arithmetic mean of the 1050 estimates. 


In tables [6^ and 6.2, by “area”, “b’dary” and “euler” we mean that in the method of 
chapter the regression has been based exclusively on one of the according data sets 
j/ 2 , 2/1 and yo (see chapter]^. Note that the case “area” corresponds to the sausage 
method from section |4.3[ By joint2 we mean the regression estimate based on both 
2/2 and 2 / 1 , and by joints the estimate based on all three data sets. 


If the set F does not belong to there is no theoretical foundation for yi and 2/0 
to be used in a regression estimate for the dimension (yet). So, strictly speaking, for 
the Koch-curve and for the Cantor dust, not only euler and joints should be marked 
“N/A” but also joint2 and b’dary. However, yi seems to admit good fits for all self¬ 
similar sets, and thus the estimates have been kept in the table. Also note that the 
Sierpihski tree is a member of yet still 2/0 cannot be used for a regression here as 
it does not represent the total variational measure Cq°‘'^{Fx-) very well. 

It seems that data given by the box-counts Ns and by the rescaled areas 2 / 2 / were the 
most reliable, as they allowed a good £t and yielded an accurate slope. The rescaled 
boundary lengths 2/ij allowed almost as good a fit, however the slope seemed to be 
systematically to low by around 0.05, especially for the square-like sets, which tended 
to have a slope too low by almost 0.15. The rescaled euler numbers 2 / 0 / were the least 
reliable data set. The goodness of its fit strongly depended on the non-arithmeticity 
of the sets, the best £t being achieved with A, the self-similar triangle. As more 
“realistic” fractals will almost certainly be non-arithmetic, there is the hope that this 
data set will perform better on non-artificial sample images. 
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images of size 3000 x 3000 



theor. 

box 

local 

area 

b’dary 

euler 

joint 2 

joint 3 


1.262 

1.254 

1.281 

1.211 

(1.274) 

(1.394) 

(1.242) 

(1.291) 


1.262 

1.270 

1.354 

1.268 

(1.237) 

N/A 

(1.252) 

N/A 



1.585 

1.564 

1.617 

1.548 

1.508 

N/A 

1.527 

N/A 

A 

A, A. 

1.585 

1.539 

1.568 

1.585 

1.555 

1.607 

1.570 

1.582 

Ik 

1.585 

1.588 

1.563 

1.583 

1.552 

1.557 

1.574 

1.564 


1.588 

1.573 

1.555 

1.585 

1.553 

1.571 

1.569 

1.569 




1.774 

1.706 

1.760 

1.732 

1.656 

1.635 

1.694 

1.715 




1.774 

1.730 

1.791 

1.744 

1.642 

1.758 

1.692 

1.670 


Ittt- 


1.774 

1.751 

1.790 

1.748 

1.667 

1.809 

1.708 

1.752 



1.794 

1.781 

1.802 

1.772 

1.710 

1.897 

1.741 

1.815 




1.893 

1.866 

1.872 

1.857 

1.723 

1.923 

1.790 

1.893 




1.893 

1.826 

1.879 

1.857 

1.715 

1.914 

1.788 

1.912 


1.893 

1.778 

1.785 

1.839 

1.605 

1.851 

1.766 

1.819 




1.896 

1.819 

1.859 

1.854 

1.690 

1.786 

1.854 

1.854 


ffi 


1.896 

1.839 

1.875 

1.863 

1.716 

1.766 

1.789 

1.945 


Table 6.1: Comparison of dimension estimates for large-sized images. 


6.3 Estimates of the Fractal Curvatures. 


Here we describe the first attempt we know of to measure fractal curvatures with 
an image analyser. Since we expect this task to be very difficult, we used only the 
highest resolution versions of the sample images (3000 x 3000 pixels). 


Like most image analysers, the algorithms from this thesis all use the pixel scale, 
i.e. the length 1 corresponds to the distance between two (horizontally or vertically) 
neighbouring pixels, meaning that on this scale the base of e.g. the Sierpihski Carpet 
has length 2980 and not 1 as assumed in Winter’s calculations. 


Since the sample sets all satisfy Sk = = s — k, hj proposition 2.2.4 we can switch 

between the two scales by multiplying by A^, where A is the scale ratio and s is the 
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images of size 250 x 250 



theor. 

box 

local 

area 

b’dary 

euler 

joint 2 

joint 3 


1.262 

1.312 

1.366 

1.130 

(1.268) 

(1.489) 

(1.187) 

(1.275) 


1.262 

1.220 

1.339 

1.228 

(1.170) 

N/A 

1.204 

N/A 



1.585 

1.571 

1.673 

1.570 

1.545 

N/A 

1.539 

N/A 

A 

A, A. 

1.585 

1.548 

1.660 

1.559 

1.519 

1.551 

1.521 

1.499 

Ik 

1.585 

1.516 

1.657 

1.553 

1.513 

1.512 

1.515 

1.505 


1.588 

1.535 

1.627 

1.545 

1.495 

1.517 

1.498 

1.537 




1.774 

1.709 

1.792 

1.721 

1.650 

1.844 

1.648 

1.749 




1.774 

1.738 

1.791 

1.733 

1.601 

1.555 

1.670 

1.556 


Ittt- 


1.774 

1.799 

1.798 

1.727 

1.616 

1.791 

1.632 

1.722 



1.794 

1.726 

1.803 

1.748 

1.704 

1.652 

1.680 

1.745 




1.893 

1.797 

1.874 

1.821 

1.735 

1.665 

1.713 

1.894 




1.893 

1.747 

1.878 

1.813 

1.610 

1.678 

1.662 

1.899 


1.893 

1.823 

1.880 

1.824 

1.607 

1.780 

1.675 

1.814 




1.896 

1.805 

1.891 

1.819 

1.739 

1.279 

1.694 

1.876 


ffi 


1.896 

1.834 

1.855 

1.825 

1.775 

1.710 

1.712 

1.895 


Table 6.2: Comparison of dimension estimates for small-sized images. 


fractal dimension: 

CliXF) = A**+*C;(F) = X'CliF). 


Table 6.3 shows a comparison between the fractal curvature measures which have been 
calculated this way and the estimated fractal curvatures, obtained by the algorithm 
from chapter]^ When comparing the fractal curvatures, one should keep in mind that 
the “unit” is pixeP, and thus curvatures are only comparable if the corresponding 
sets have the same dimension. Note also that one should make sure that each set has 
equal base length, as is the case here (all bases have length 2920 pixels). 


- J 

A typical relative error lies at around 10%, however, for the rectangular-like sets, 
seems to be systematically overestimated. This corresponds with the observation from 
the previous section, where the measured growth of the rescaled boundary lengths 
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dim 

'^0 

ci 

ci 



1.585 

? 

0.0 

? 

145556 

? 

631342 

A 

AAAA 

1.585 

-13197 

-11109 

117230 

124064 

564100 

568985 



1.585 

? 

-9214 

? 

106250 

? 

498975 


1.588 

-9843 

-7733 

100416 

94940 

487649 

448388 


n 


1.893 

-58716 

-59997 

262770 

386937 

4900200 

4861736 




1.893 

-50742 

-48113 

260960 

339448 

4871275 

4627544 


i: 


1.896 

-54228 

-42878 

242166 

292915 

4644054 

4306752 


i 


1.896 

-60811 

-65169 

271565 

371606 

5207828 

5843391 




1.774 

? 

-30530 

? 

280567 

? 

2056714 




1.774 

? 

-6699 

? 

278438 

? 

2093975 


□ o 

13 o 


1.774 

? 

-9380 

? 

351984 

? 

2754060 


Table 6.3: Estimated Fractal Curvatures. In each row, the upper number is the 
theoretically calculated curvature, whereas the lower number is the estimate returned 
by the algorithm from chapter A question mark denotes curvatures that have not 
been calculated yet, and an “N/A” mark means that for the corresponding fractal 
either Cq (F) = 0 or the asymptotic of yoj is not good enough. 


yij was too low, also especially for the rectangular-like sets: Recall the calculation 


:= ^^ - tj_i) 


tm to 




of the k-th fractal curvature estimate on page 52 If s is assumed higher than a least- 
squares £t of yij vs. Xj would suggest, the convexity of the exponential function will 
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yield a relatively too high value, since the Xj are all negative. The estimator exp(Zi)i) 
shows similar behaviour. 


The last four rows in table [6^ show the two pairs “Window and Gate” and “Ml and 
M2”. In both pairs, the two sets have the same dimension, but the first one seems 
more lacunar. This is reflected in a lower Minkowski-content, as has been already 
calculated for the hrst pair. The author’s conhdent guess is that it is also the case 
for the second pair. In any case, lacunarity seems to have the effect of reducing all 
of the curvature measure estimates. As this behaviour is parallel to the dimension 
estimates, the absolute values of the curvatures might not be suitable to discern 
different fractals of equal dimension. However, their size relative to each other yields 
an important geometric invariant (see nxt section). 


6.4 Characterization of Q^-sets beyond dimension 


Specific fractal curvatures. For a systematic categorization of fractal sets, ge¬ 
ometrically invariant characteristics are useful. Here, by a characteristic we simply 
mean a functional 

F 

on the class IK'’* of compact sets in M'’*. F is called 

• motion invariant if F{g{K)) = F{K) for all K G IK'’* and for all euclidean 
motions g on M'’*, 

• homogeneous of degree k if F{\K) = X^F{K) for all A G (0, oo), and 

• scaling invariant if F is homogeneous of degree 0. 


For example, all definitions of fractal dimension satisfy motion and scaling invariance. 
We now define the following characteristics: 


Definition 6.4.1 Let A G IK'’* be such that cl{A) exists for all k G {0,..., d}, and 

- J 

assume C^{A) > 0. Then for each fc G {0,..., d — 1} 




ci(A) 

0^(A) 


is called the k-th specific fractal curvature of A. 
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Figure 6.15: Local observations of the Sierpinski gasket (left) and A (right). 

Note that C^{F) is always positive for self-similar sets satisfying the OSC as shown 
by Gatzouras in |Gatnn] . so is well-dehned for all Q^sets. 

Since the fractal curvatures are motion-invariant and homogeneous for all fc G {0,..., d— 
1}, so are the specihc fractal curvatures. There are the following two important cases: 


Proposition 6.4.2 Assume that ^^(A) is well-defined, and let A > 0. 

1. If A E TZ^, then 

EkiXA) = X'^-^EkiA). 

2. If Sk{A) = s — k for all k E {0,..., d}, then 

2fc(Ay4) = Ek{A). 

Proof: This is a simple consequence of the Sk -l- k- homogeneity of Cl, the Sk being 
0 for sets in the convex ring. □ 

Local observations. Let be the usual geometrical similarity relation, re¬ 
stricted to Q'^. If for two sets F,G E Q'^ their binary representations F and G 
are given and one has to decide whether or not F G, then F and G can be rescaled 
to have the same diameter, and their (average) Minkowski-contents can be estimated. 
If they differ signihcantly, this could be evidence against F ~ G. 

However, in many applications F and G will only be local representations, i.e. binary 
images of F fl Wi and G fl W 2 , respectively, where the Wi are (e.g. rectangular) 
bounded observation windows. If F is to be compared to G, one faces the problem 
that along with the position of the window also the amount of white space will vary. 
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and thus lacunarity analysis by the box-gliding algorithm (chapter and any total 
curvature C^{F) will not be a reliable source of information. (Note that it might 
happen that F n Wi does not have polyconvex parallel sets anymore even though F 
does, and that thus C^(F fl Wi) might not be dehned any more; but this problem 
gets lost in the discretization procedure anyways.) However, the (average) Minkowski- 
content 6 * 2 ( 7 ^) can serve as a normalization factor for the other fractal curvatures, 
i.e. the specihc fractal curvatures can still be calculated. If they are signihcantly 
different then F and G might still be distinguished from each other even though the 
dimension appears to be the same. 


6.5 Discussion 


Dimension estimates. The tests on the sample images suggest that the overall 
accuracy of dimension estimates via the simultaneous regression on all data sets “euler 
number”, “boundary length” and “area” is comparable to the accuracy of the sausage 
method and the box-counting method. More data do not lead to more accuracy here: 


The measured growth of the boundary length of the parallel sets as £ 0 

is slightly lower than the theory would suggest (see column “b’dary” of ta¬ 


ble 6.2). Measurements based on the algorithms of references |KSS06j and 
|OMOOj both qualitatively yielded almost the same results; the latter algorithm 
showed slightly slower growth. A partial reason for this might be the convexity 
of the plot Hi = log against x = — log£, which does not die off quickly 

enough as x approaches inhnity. Cutting off the data before the break condition 
Cq°'‘^{F^) < 2 applies might help getting rid of some of the negative bias that 
the slope of the plot has, but the overall stability of the estimate will become 


worse. 


• For the measurec02A^(Fe) — Co(F7) to be close to the total variational measure 
Cq“”(F£), the holes of F^ need to be convex (also see appendix). But even if 
they are convex, they might exhibit stairlike behaviour due to arithmeticity, 
which does not allow for a good regression £t. Of the considered sets, only the 
self-similar triangle allowed for a good fit of all three data sets. 

^Recall that N{FA is the number of connected components of 
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Curvature estimates. The accuracy of fractal curvature estimates strongly de¬ 
pends on the accuracy of the dimension estimate. As noted before, the data set of 
Cq'^^(Fs) can seldomly be measured properly, and thus it should be included in the re¬ 
gression analysis only in special cases. Note however that even if Cq‘^^(Fs) is excluded 
from the regression estimate of the dimension, then it is still possible to calculate the 


0-th fractal curvature estimate ro(F) (compare equation (5.12)). 

The functional dependence of the curvature estimates on the dimension estimate can 


be read off the following alternative formulation of equation (5.12): 

1 


p(m) _ 

^ k — 


tr. 


to 


sgn(Cfc(Fg-.,)) exp {ykj 


SXi 


(tj 


ti-l) 




If the slope s is higher or lower than the plot of i/kj against xj suggests, there will 
be higher and lower arguments to the exponential function, and since it is convex, 
positive deviations will come out stronger than negative deviations. Thus a wrong 
dimension estimate will always result in too high values for no matter if the 

dimension estimate itself will be too high or too low. As this behaviour is the same 
for all k G {0,1,2}, one might hope that the negative influence of a bad dimension 
estimate will cancel out for the specific curvatures So(A) and Si(F), and that thus 
So(F) and Si(F) are less susceptible to wrong dimension estimates. 


Specific curvatures. As pointed out before, it is hard to compare two sets by 
their (average) Minkowski-content, if one cannot be sure that the same scale should 
be used on both sets; in this case, the specific curvatures seem more appropriate, 
though the Minkowski-contents should be kept in mind. 


Table 6.4 shows the specihc curvatures of the sample sets in the right columns, where 
in each row the upper value is theoretical and the lower value is estimated from 
the binary image. The general observation is that the specific curvatures increase 
as the dimension of a self-similar set decreases from 2 to 1. Also note that for the 
observed sets, if two fractals have indistinguishable dimensions then their specific 
curvatures may be different; however, in most cases like this the specific curvatures 
will be too close together for an image analyser to depict the difference with high 
enough confidence. There are cases, however, where the specific curvature estimates 
can distinguish sets where the dimension estimates cannot: As an example, compare 


the Tripet to the self-similar square □ in table |6.4[ The box-counting estimates are 
almost the same, but Si is rather different. 
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dims 

dim 




1.254 

1.262 

N/A 

N/A 

(0.412) 


1.270 

1.262 

N/A 

N/A 

(0.373) 

A 

1.539 

1.585 

-0.0234 

-0.0195 

0.208 

0.218 



1.564 

1.585 

0 

0 

? 

0.230 



1.573 

1.588 

-0.0202 

-0.0172 

0.206 

0.211 



1.588 

1.585 

? 

-0.0185 

? 

0.213 


Trrr- 

XEZ 


1.706 

1.774 

? 

-0.0148 

? 

0.136 


P 


1.730 

1.774 

? 

-0.00319 

? 

0.1329 


EH m 

EH eh 


1.730 

1.774 

? 

-0.00340 

? 

0.1278 



1.778 

1.893 

? 

-0.00715 

? 

0.0675 




1.781 

1.794 

? 

-0.00577 

? 

0.1114 




1.819 

1.896 

-0.0117 

-0.0100 

0.0521 

0.0680 




1.826 

1.893 

-0.01041 

-0.01039 

0.0536 

0.0733 


m 


1.839 

1.896 

-0.0117 

-0.0112 

0.0521 

0.0636 

n 

1.866 

1.893 

-0.01183 

-0.01234 

0.0536 

0.07958 


Table 6.4: Estimated specific fractal curvatures: The first two columns contain the 
estimated dimension via box-counting and the theoretical value, respectively. In the 
right two columns, the upper value of each row is theoretical and the lower value is 
the returned estimate. The Cantor Dust and the Koch Curve are marked “N/A” 
since they are not in 










Concluding Remarks 


In most of the literature on estimates of fractal dimension, the slope of a logarithmic 
plot of Cd{Fj.) against r earned almost all attention, whereas the intersection point of 
the htted line with the vertical axis was mostly ignored. Furthermore, in the context 
of fractal analysis the author has found almost no occurence of any other generalized 
volumes than Vd (= Cd) in the literature; The only exception was B. Mandelbrot’s 
exposition |Man94j on “gap-lacunarity” of Fractals on the one-dimensional unit in¬ 
terval, which is essentially the 0-th fractal curvature Cq. This thesis shows that, at 
least for the special case of Q'’*-sets, a large amount of information is ignored if only 
the d-dimensional volume of the parallel sets is analysed. 


The theory of fractal curvatures measures is still at its beginning stage, and only for 
self-similar sets there are numerical expressions for their calculations. Their dehnition 
can be somewhat extended, namely to sets whose parallel sets are hnite unions of sets 
of positive reach; However in that case little to nothing is known about the scaling 
exponents Sk, or even if the growth of Ck{Fj.) can be compared to for some s. But 
there does seem to be a similar behaviour of self-Oj^ne sets in terms of the scaling 


behaviour of the curvatures of parallel sets, as Figure 6.16 shows. 


An example for the case where Ck{Fr) ~ F does not hold true for any s G M is 
a brownian path i?t([0,l]) on the time interval [0,1] in For this statistically 
self-similar set. 


EC'2(Rt([0,l]),) 
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log r I ’ 


EC'i(Hi([0,l]), 


TT 


(r;0), (6.1) 


2r log r 

the second asymptotic having been proved only recently in |RSSn5] . Numerical sim¬ 
ulations of EC'o(i?i([0, l])r can be found in [RSMOTj . as well as an empirical conhr- 
mation of equations ( 6.1[ ). In this case, regression analysis on the log-log scale does 
not seem apropriate any more, and other means have to be found to measure Fractal 
Curvature. 
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(a) Only two affine maps on generate this attractor. 



(b) The usual plot of ykj against Xj, k G {0,1, 2}. The fit seems to be just as good as for an Q'^-set. 

Figure 6.16: The Dragon, a self-affine set 
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Appendix 


Estimating Cq^'^^Fs) 

In general, with an image analyser only the Euler-number Co{F^) of the parallel 
sets can be determined; for the regression analysis of chapter]^ however, 
is needed in order to hnd the scaling exponent Sq. Since in 2D, the Euler-number 
equals the number N of connected components minus the number Q of holes, one 
could think that the total variational measure Cq°'^{F^) would equal N + Q, but this 
is false in general. The following lemma gives a sufficient condition for this to be true. 


Lemma 1 Let F G TZ^, and let N be the number of connected components of F and 
Q the number of holes of F (i.e. the number of connected components of F'^ minus 
1). Moreover, assume that the closures of the components of F^ are disjoint. Then 


1 . |C'o(F)| < A + g<C'o™"(F). 

2. Let bd F = IJ^i be the decomposition of the boundary of F into its connected 
components Tj. If 

= 1 for all i, (6.2) 

then 

C-ar^F) =N + Q. 


Proof: Since C'o(F) = N — Q the left inequality in 1 is clear. Now we use the result 
M = N + Q, where M is the number of connected components of the boundary of F\ 

To prove this, start off by labelling the components of F and of F'^ in the follow¬ 
ing way: Assign 0 to the only unbounded component. Assign 1 to all remaining 
components of F that have a boundary with 0, and proceed with assigning i to the 
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remaining components that have a boundary with the components labelled i — 1 until 
all components of both F and have been labelled. By assumption, no component 
of has a boundary with any other component of and thus components of F'^ 
have even labels, whereas components of F will have odd labels. There will be a 
maximum label m, since F G Now for all odd k define F^ to be the union of all 
components with odd labels up to k. Then the number of components of the bound¬ 
ary of Ffc will be the equal to the number of components with label < (fc -|- 1) minus 
1, as can be checked by induction on k. In particular, the boundary of F = Fm-i 
will have as many components as the total number of components of both F and F'^ 
together, minus 1. But this is the number of connected components of F plus the 
number of holes of F. 


Since always FQ“^(rj) > 1, the right inequality in 1 follows, and 2 is just a consequence 
oi M = N + Q. □ 


Note that the condition C'Q“'’(rj) = 1 is equivalent to the interior of the curve Tj 
being convex. Thus out of all sets considered in section 6.1, equality (6.2) applies in 
the following way to our sample sets: 


• It does apply to Sierpihski Gasket and Carpet, Tripet, A, Window and Gate. 


It does not apply to the parallel sets of the Cantor dust and Koch curve, since 
the parallel sets are not polyconvex. 

For the set Ml, in fact Go“'’(rj) = 3, so, for small £, the value Cq°‘^{F^) will be 
3 times as high a.s N + Q. This will reflect in a much too low estimate of the 
0-th fractal total variational curvature, but it should not affect the estimate of 
c;(Mi). 


For the modified Sierpihski Carpet, assuming equality (6.2) is only a little mis¬ 
take: Here 

0 < CriF,) - (iv + 0) < J 

for all £ > 0, as only the above opening is causing trouble. 


• Unfortunately, for the set □ the value N + Q is not close to CQ“'’(Fe) at all: 
Depending on e, the holes will seldom be convex and have the total varia¬ 
tional curvature 1 as demanded by the above condition. Instead, often they 
are L-shaped, or cross-shaped, thus having total variational curvature 1.5 or 3 
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respectively, and other shapes are possible. Thus, for this set useEuler is best 
set to false. 

• The Sierpihski Tree is another example where the 0-th total variational curva¬ 
ture cannot be measured; all parallel sets have one connected component and 
no holes, whereas the boundary is obviously very curved. 


Calculation of the total Fractal Curvature of A 


We need the curvature scaling functions 


N 


Rt{E) = ft(A,) - l(0«|(E)Ct((S.A), 


(6.3) 


2=1 


for k G {0,1,2} on the interval (0,1]. It turns out that the Rk{s) are piecewise 
polynomials of degree k whose transition points are the discontinuities of the indicator 


functions at ri = H, r2 = tt ’"a = if’ radius of the in-circle u (see 


20 


16 


hgure 6.10). The inner triangle is in fact congruent to the upper (red) smaller copy, 
and thus we can compute f |Wik07j ]) 


oj = 2r2 


\ac 


a-\-h + c 


4 

41' 


Now we determine the curvature scaling functions directly via formula (6.3). Note 
that for 0 < £ < ra the inclusion-exclusion formula says that 


Rkie) = 


|/|>2 Vie/ / 


the sum being taken over all subsets / C {1,2,3} with at least two elements, and 
that for 0 < £ < ca the intersection 


(^i(A)),n(^ 2 (A)),n (53(A)), 


is empty. 

SO that in 

this case R 

k{e) = - 

Ei^j Ck (Eij) where 

' Eij — 

5i(A),nA(A),. 

We hnd 







Cf){Ei2) 

= 1, 


A (As) 

= 1, 

A (As) 

= 1, 

C'i(Ai2) 

= (1 + 

76+7 \ 

2 

A (As) 

= (2 + A^)e, 

A (As) 

= (S + A')^, 

C2{Ei2) 

= (1 + 

I + 7) A 

^(As) 

= (2 + f + A^V 

A (As) 

= (3 ^ -f- Cl)£^ 
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Now, for each k G {0,1, 2} we use the representation 


Rk{e) = < 




ri 

< 


< 

1 

Ck{A,) 

-Afc((5iA),), 

r2 

< 

£ 

< 

ri 

Ck{A,) 

-Afc((5iA),)-Afc((^2A),), 

rs 

< 

£ 

< 

r2 

Ck{A,) 

-Yl=iCk{{SkA\), 

00 

< 

£ 

< 

r3 


Ck{Eij), 

0 

< 

£ 

< 

00 


and putting 


A-= — = 0.24, 
2 


S :=a + 6 + c = 2.4 
for the area and the circumference of the big triangle we arrive at 


Rnis) — ■* 


^ § + 7ie, 


(1 - ’’Of, 


-Ri(^) — ■* 


1, 

ri < £ < 1 

0, 

r 2 < e < ri 

-1, 

rs < e <r 2 

-2, 

oo < e < Ts 

-3, 

0 < e < 00 

■Of- 

Tie, ' 


ri < e < 1 
r 2 < e < ri 
rs < e <r 2 


0 < e < uj 


R2ie) = < 


(1 - ri 

(1 - ri - r 2 - rs)! - 27re, oo < e < 
-(6 + 27r)e, 

y4 + S'e + 7r£^, 

(1 - r^A + (1 - ri)S'£, 

(1 — — r|)y4 + (1 — ri — r2)Se — 7r£^, 


ri < £ < 1 

r 2 < e < ri 
rs < e <r 2 


(1 


r‘^)A + (1 — ri — r 2 — r 3 )S — 27r£^, OJ < e < 


■(6 + 27r)£^, 


Now the fractal curvatures C'^(A) can be computed as in theorem 
author used a shortcut-formula from |Win05j . 


0 < e < 00 

but the 


2.3.7 
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Signed measures 


Definition 2 (signed measure) Let {X, S) be a measurable spaee. If the mapping 

/i : S —>■ K. 


satisfies 

( oo \ oo 

i=l ) i=l 

for a sequence of disjoint sets Ai G S, then it is called a signed measure. 

Definition 3 (positive, negative and total variational measure) Let n be a signed 
measure on the measure space (X, S). Then the positive variational measure, negative 
variational measure and the total variational measure are defined respectively by 

ia'^{A) := sup{/i(-B) : B G E, B C A} 

/i" (A) := - inf{/i(5) :B eE,BC A} 

/i"“"(7l) := fi+{A)+fi-{A). 

which are non-negative measures on (X, S). 

By the Hahn decomposition theorem, fjA — p,~ = fi for every signed measure p,. 
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